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ABSTRACT 

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz muhi-band 
imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 
4-meter telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point- 
spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure 
variation of this position-dependent PSF change is significant. We overcome this challenge by modeling 
the PSF separately for individual exposures and CCDs with principal component analysis (PCA). 
We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high 
fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing 
effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates 
sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic 
atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear 
calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology- 
dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this 
varying covariance is critical in our parameter estimation. Our current non-tomographic analysis 
alone constrains the JIm — cts likelihood contour tightly, providing a joint constraint of f2j\/ — 0.262 ± 
0.051 and erg = 0.868 ± 0.071. We expect that a future DLS weak-lensing tomographic study will 
further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear 
more efficiently breaks the VLm — erg degeneracy. Combining the current results with the Wilkinson 
Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain VIm = 0.278 ± 0.018 and 
0-8 = 0.815 ±0.020. 

Subject headings: cosmological parameters — gravitational lensing: weak — dark matter — cosmology: 
observations — large-scale structure of Universe 



1. INTRODUCTION 

Weak gravitational lensing from large-scale structures 
in the universe, often called cosmic shear, allows one 
to address a number of critical issues in modern cos- 
mology. Its application encompasses the study of the 
universe's matter density and its fluctuation, probes of 
the footprints of non-Gaussianity in the primordial den- 
sity fluctuation, constraints on dark energy and its evo- 
lution, tests for modified gravity, etc. The consensus 
on the critical role of cosmic shear studies triggered 
quite a few optical surveys such as the Canada-France- 
Hawaii- Telescope Legacy Survey (CFHT-LS; Hoekstra 
et al. 2006, Semboloni et al. 2006; Fu et al. 2008), 
the Red-sequence Cluster Survey (RCS; Gladders et 
al. 2003), the Cerro Tololo Inter- American Observa- 
tory (CTIO) Lensing Survey (Jarvis et al. 2006), the 
Garching-Bonn Deep Survey (GaBoDS, Hetterscheidt et 
al. 2007), the VIRMOS-DESCART survey (VIRMOS, 
Van Waerbeke et al. 2005), the Deep Lens Survey (DLS, 
Tyson et al. 2001, Wittman et al. 2006), etc. The 
current surveys include the Dark Energy Survey (DES, 
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The Dark Energy Survey Collaboration 2005), the Kilo- 
Degree Survey (KIDS; Verdoes Kleijn et al. 2011), the 
Panoramic Survey Telesope and Rapid Response System 
(Pan-STARRS, Kaiser et al. 2010), etc. Next-generation 
weak-lensing projects are the Euclid mission (Laureijs 
et al. 2010), the Wide Field Infrared Survey Telescope 
(WFIRST; Green et al. 2011), and the Large Synoptic 
Survey Telescope (LSST, LSST Science Collaborations 
et al. 2009). 

Needless to say, great effort should be given to control 
of systematics in both shear and photometric redshift 
measurements for these future surveys. The unprece- 
dentedly small statistical errors will bring revolutionary 
advances to cosmology only if progress in shear cali- 
bration and control of catastrophic errors in photomet- 
ric redshift estimation parallels the increase in statisti- 
cal power. The recent shear estimation challenges such 
as the Shear TEsting Programme (STEP, Massey et al. 
2007; Heymans et al. 2006), the GRavitational lEnsing 
Accuracy Testing (GREAT, Bridle et al. 2009; Kitching 
et al. 2012), etc. are concerted efforts to quantify bias in 
the current popular shear estimation methods and also 
to identify the limitation of the current weak-lensing sim- 
ulation methods. Similar efforts toward improvement of 
photometric redshift estimation, albeit less mature, are 
also underway (e.g., Hildebrandt et al. 2010). 

Both sky coverage and depth must be carefully bal- 
anced to maximize the scientific return from future cos- 
mic shear surveys. Large sky coverage is needed to min- 
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Figure 1. Survey area and depth of various optical surveys. The 
red Une represents the AQt = constant locus, where A, fi, and 
t are the primary mirror area, field of view, and exposure time, 
respectively. DLS is the deepest optical survey to date among the 
current > 10 sq. degree surveys. Depth is compared either in the 
B ovi band. 

imize the contribution to the error from the sample vari- 
ance. Deep imaging is required to detect and measure 
the shapes of high redshift sources, which ahows us to 
probe the evolution of the cosmic structure over a signif- 
icant fraction of the age of the universe. The Deep Lens 
Survey (DLS; Tyson et al. 2001, Wittman et al. 2006) 
is designed as a precursor to these next generation cos- 
mic shear surveys with emphasis on the latter, reaching 
a mean source redshift of 2 ^ 1 over 20 sq. degrees us- 
ing the two National Optical Astronomical Observatory 
(NOAO) Mayall and Blanco 4-meter telescopes. Fig- 
ure [1] shows the comparison of DLS sky coverage and 
depth with those of other optical surveys. DLS is the 
deepest optical survey to date among the current ^10 
sq. degree surveys. Galaxy populations are dominated 
by faint blue galaxies when a survey reaches or exceeds 
the depth of the DLS. As no ground-based cosmic shear 
study with a comparable depth has been presented, the 
current cosmic shear analysis with the DLS is an impor- 
tant experiment, testing whether the shapes of the faint 
blue galaxy population smeared by atmospheric seeing 
can be reliably used for cosmic shear. We augment this 
experiment using image simulations with real galaxy im- 
ages. In addition, our seeing-matched photometry from 
the deep BVRz imaging in conjunction with a large spec- 
troscopic sample allows us to stabilize our photometric 
redshift estimation and to identify where potential sys- 
tematic errors lie in our results. Reliable photometric 
redshifts are pivotal not only in the interpretation of the 
cosmic shear signal, but also in future application of the 
measurements to weak lensing tomography. 

This paper is the first in a series of our DLS cosmic 
shear publications. Here we mainly focus on the DLS 
systematics induced by the point-spread- function (PSF), 
the removal of the systematics with our principal com- 
ponent analysis (PCA) and "StackFit" methods, and the 
two-dimensional (non-tomographic) analysis of the DLS 



cosmic shear signal. 

The structure of this paper is as follows. In fj2l we de- 
scribe our DLS data and analysis method including our 
detailed PSF modeling and shear calibration efforts. The 
theoretical background of cosmic shear and our system- 
atics control is presented in fJ31 We discuss the study of 
cosmological parameter constraints in 21 ^^^d conclude in 

m 

2. OBSERVATIONS 
2.1. Data 

The detailed description of the DL^ can be found in 
Wittman et al. (2006; 2012). Below we provide a brief 
summary of the survey and its data. 

The DLS covers five 2° x 2° fields (hereafter F1-F5). 
Fl and F2 are in the northern sky, and observed with the 
Kitt Peak Mayall 4-m telescope/Mosaic Prime-Focus Im- 
ager (MuUer et al. 1998). F3, F4, and F5, which are in 
the southern sky, were observed with the Cerro Tololo 
Blanco 4-m telescope/Mosaic Prime-Focus Imager. Ta- 
ble 1 lists the coordinates of the five fields. 

Each Mosaic Imager provides a ~ 35' x 35' field of view 
with a 4 X 2 array of 2 kx 4 k CCDs (~0'.'26 per pixel). 
We divide each 2° x 2° DLS field into a grid of 3 x 3 
array. Each 40' x 40' subfield, slightly larger than the 
camera field of view, was covered with dithers of -200". 
The DLS data consists of 120 nights of S, V , R, and 
z imaging. A priority was given to the R filter, where 
we measure our lensing signal, whenever the seeing was 
better than ~ 0'.'9. The mean cumulative exposure time 
in R is about 18,000s per field whereas it is about 12,000s 
per field for each of the rest of the filters. The typical 
exposure time per visit is about 900s. 

2.2. Reduction 

We applied initial bias, flat, and geometric distor- 
tion correction to the DLS data with the IRAF package 
MSCRED. External astrometric calibration was performed 
by matching astronomical objects in each exposure to the 
USNO-Bl star catalog using the msccmatch task. The 
residual uncertainty in the global coordinate system rel- 
ative to the USNO-Bl catalog is less than 0.01". The 
mean rms error per object is ~0'.'3. The limiting factor 
for this scatter per object is believed to be the internal 
accuracy of the USNO-Bl catalog. Internal astrometric 
calibration between different epoch data was carried out 
using the common high S/N stars present in the overlap- 
ping region. Precise registration is essential in precision 
weak-lensing analysis because a small -0.5 pixel error 
can create a noticeable correlation of object ellipticity 
over a large scale. We verify that the mean rms error per 
object is less than ~0.1 pixel and the scatter is isotropic, 
which indicates that the scatter is dominated by photon 
noise. 

We found an initial non-negligible (10—20%) residual 
fiatfielding error in the final stack image after the appli- 
cation of the sky-fiat correction. This is further refined 
to the 2—5% level using the "iibercal" method (Padman- 
abhan et al. 2008). Interested readers are referred to 
Wittman et al. (2012) for details of our "ubercal" imple- 
mentation and performance. 

* http:/ /dls. physics. ucdavis.edu. 
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Table 1 

DLS Fields and Data 



Field Name RA DEC Median Seeing (R) 

^source "^source 

(") (per sq. arcmin) 



Fl 00:53:25 +12:33:55 

F2 09:18:00 +30:00:00 

F3 05:20:00 -49:00:00 

F4 10:52:00 -05:00:00 

F5 13:55:00 -10:00:00 




2000 4000 6000 8000 

X (pixels) 



Figure 2. Example of spatial variation of DLS PSF. Although 
this particular pattern is observed on 24 February 2001 from the 
Blanco telescope, a similar degree of PSF variation complexity is 
commonly present in all of our DLS data. Each "whisker" shows the 
direction and magnitude of the stellar ellipticity at the location by 
its orientation and length, respectively. The red stick in the middle 
shows the size of 10% ellipticity [i.e., (a - b)/{a + b) = 0.1]. The 
eight shaded rectangles depict the eight CCDs of the camera. Here 
we did not clean up outliers (e.g., cosmic-ray hit stars, binary stars, 
etc.), and they do not represent real PSFs. 

Our team has developed two pipelines (Pipeline I and 
II) for the creation of the final mosaic. Pipeline I is 
optimized for photometry and consists of independently 
implemented standalone programs (Wittman et al. in 
preparation). It performs PSF-matched photometry to 
minimize the systematics in photometric redshift estima- 
tion (Schmidt et al. in preparation). Pipeline II is opti- 
mized for weak lensing and controls the flow of the SCAMP 
and SWARP program^ We process only R-band data with 
this second pipeline. These two pipelines share the above 
procedures, but differ in that the weak-lensing pipeline 
uses the subset (with better seeing and less astrometric 
issues) of the DLS data and creates a large 2° x 2° mosaic 
image per field whereas Pipeline I produces nine (3 x 3) 
subfield images to cover each 2° x 2° field. In i i2.3l we de- 
scribe this weak-lensing pipeline in detail in the context 
of the PSF reconstruction. 

2.3. PSF Reconstruction 
^ available at |http: / /www. astromatic.net | 
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The spatial variation of the PSF is substantial and 
complicated for both the Mayall and Blanco telescopes. 
An example of this PSF pattern is displayed in Figure [2j 
Although this particular pattern is observed on 24 Febru- 
ary 2001 from the Blanco telescope, a similar degree of 
PSF variation complexity is commonly present in all of 
our DLS data. It is difficult to interpolate the variation 
over the entire focal plane with a single set of polynomi- 
als. Thus, polynomial interpolation should be limited to 
a smaller area, where the variation is slow and tractable. 
Hence, we choose to model the PSF variation on a CCD- 
by-CCD basis. This chipwise approach was investigated 
by Jee et al. (2011) for the LSST, where the smaU /- 
ratio of the optics makes the potential aberration highly 
sensitive to CCD flatness, giving rise to a sudden, no- 
ticeable jump in PSF patterns across CCD boundaries. 
For the Mayall and Blanco telescopes, we often found 
a somewhat smaller, but clear discontinuity across the 
CCD gaps, although in principle the relatively large /- 
ratio of the two telescopes should make the CCD-to-CCD 
flatness much less important. 

As most lensing signals come from distant, faint galax- 
ies, which sometimes are not even detected in single ex- 
posures, these source galaxies are commonly examined 
after multi-epoch data are combined to produce the deep 
stack image (i.e., single 900s exposure vs. cumulative 
18,000s exposure). Therefore, it is important that the 
PSF modeling closely mimics the image stacking proce- 
dure (e.g., offsets, rotations, geometric distortion correc- 
tions, etc.). Figure [3] schematically illustrates how image 
stacking complicates the PSF pattern. After stacking is 
performed, across the image boundaries of input frames 
we often observe a discontinuous change of PSF as dis- 
played. This discontinuity prohibits us from interpolat- 
ing PSFs based on the information obtained only from 
the final stacked image. Hence, in our DLS weak lensing 
analysis, the PSF modeling is performed with the follow- 
ing two steps to address the issue. First, we construct a 
PSF model for each CCD image using a PCA method. 
Then, the PSF on the final mosaic image is computed 
by weighted combination of all contributing PSFs from 
each CCD image. Below we provide the details for each 
step. 

2.3.1. Step 1: PSF Modeling with. PCA for each CCD image 

A 2° X 2° mosaic image for each DLS field is created 
with the SCAMP/SWARP software. The SCAMP program au- 
tomatically refines WCS headers of images by first cross- 
identifying astronomical objects with external standard 
catalogs and then by tweaking the WCS information of 
each header in such a way that internal consistency is 
maximized. Because the astrometric solution is already 
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Figure 3. Complication in PSF modeling due to image stacking. Weak-lensing analysis is typically performed on a stacked image, which 
often exhibits sharp PSF discontinuities. The figure schematically shows how this complication arises. When we combine the two images 
in the left panel, the resulting PSF (right) possesses abrupt ellipticity changes across the boundaries of input frames. 

obtained in the photometric pipehne to the weak-lensing 
precision, we feed the Pipehne I catalogs into SCAMP as 
an external catalog. 

The Swarp program utilizes the series of these refined 
WCS headers to define the global WCS for the final mo- 
saic. Then, the input images are resampled and com- 
bined to create the final mosaic. We use the LanczosS 
interpolation kernel, which mimics the ideal sine kernel 
and is known to suppress the correlation between pixels. 
We estimate that the covariance between adjacent pixels 
is about 7% of the variance. This inter-pixel correlation 
leads to underestimation of both photometric errors and 
shape errors. The slight shift in shape errors also changes 
the weight in our shear correlation computation. In prin- 
ciple, we can remedy the situation by increasing our rms 
map to compensate for this underestimation. However, 
we conclude that this step is unnecessary because the 
resulting change in weight distribution is small and well 
within the interval of the shear calibration marginaliza- 
tion (g2X2|). 

What we should potentially be concerned about is the 
systematics (multiplicative) in shear calibration. The 
inter-pixel correlation somewhat smears the galaxy pro- 
file and on average circularizes the shapes. Fortunately, 
since we use the same LanczosS kernel in image simu- 
lations for our shear calibration f H2.5.2[) . the resulting 
multiplicative factor already includes this effect. 

The Swarp program provides an option to keep the in- 
termediate resampled images (hereafter RESAMP images) . 
We use these RESAMP images to identify stars and model 
PSFs because they are properly rotated, shifted, and 
distortion-corrected. For some unknown reason, some 
frames are found to possess rather large (> 0.2 pix) sys- 
tematic offsets with respect to the stacked image. In 
addition, the PSF of some frames are significantly larger 
than our criterion (FWIIM=1"). About 5% of the data 
fall into this group, and we exclude these images for the 
creation of the final stack. 

The 2° X 2° mosaic image for each field consists of more 
than ~ 1,200 CCD images. Consequently, we need an 
automatic procedure to select high S /N isolated stars and 
apply PCA to them. Our star-selection algorithm relies 
on the size versus magnitude relation with some impor- 
tant fail-safe procedures. The algorithm starts with an 
initial guess of the half-light radius and magnitude range 
of the "good" stars. Of course, because of the varia- 
tion in telescope seeing and exposure time, the stellar 
locus shifts exposure by exposure. Thus, we search the 



Figure 4. Variance vs. number of PSF basis functions 
(eigenPSFs). To determine the number of basis functions for a 
compact description of the PSF, we examine fractional data vari- 
ance for different number of basis functions. For H ST / ACS PSFs, 
we observe that the growth slows down notably after ~20 (Jee et 
al. 2007). For the PSF of the 4-m Mayall/Blanco telescopes, this 
happens at ~5. The simpler profile of the ground-based PSF (as op- 
posed to complex, diffraction limited PSF of HST) requires fewer 
basis functions. However, because of the larger FWHM variation 
(i.e., atmospheric seeing), the total variance remains slightly lower 
than in the case of HST /ACS (-96% vs. -99% at 20 ). The PSF 
reconstruction of the Mayall/Blanco Telescopes does not show any 
significant difference in quality as long as the number of eigenPSFs 
is ^ 5. In the current study, we choose to keep 20 eigenPSFs. 



two-parameter space iteratively for the stellar locus in 
the half-light radius range from 1.4 pixels to 5.5 pixels 
and the magnitude range whose minimum value (maxi- 
mum flux) is adjusted depending on the saturation level 
of the input frame. We discard stars whose SExtractor 
(Bertin & Arnouts 1996) flags are not zero. The resulting 
"clean" stars are used to derive the principal components 
(eigenPSF), and the coefficients (i.e., amplitude along 
the eigenPSF) are computed. To determine the number 
of basis functions for a compact description of the PSF, 
we examine fractional data variance for different number 
of basis functions (FigureH]). The total variance does not 
increase rapidly after five, and thus the choice is some- 
what arbitrary. We choose to keep 20 eigenPSFs, which 
accounts for -96% of the total variance. 
After we obtain these 20 eigenPSFs, the fc*'' star image 
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is decomposed as 

Ck{i,j)= X! ^k7iPn{i,j)+T{i,j), (1) 

where Ck{i,j) is the normahzed pixel value of the /c*'' 
star image at the pixel coordinate Pn is the n}^ 

eigenPSF, akn is the projection of the k*^ star in P„, 
and T is the mean PSF. Because P„'s are orthogonal 
to one another, one can determine atn by multiplying 
the corresponding eigenPSF to the mean-subtracted star 
image. 

Approximately 50-200 stars are available per CCD per 
exposure depending on galactic latitude, and we fit 3rd 
order polynomials to the spatial variation of the coeffi- 
cients to enable interpolation at any arbitrary position 
within the CCD. When we experiment with 4th order 
polynomials instead, the interpolation becomes occasion- 
ally unstable for some frames, where the number of high 
S/N stars is not sufficient. In addition, we find that the 
interpolation by 2nd order polynomials slightly underfit 
the spatial variation with respect to the 3rd order poly- 
nomial result, increasing the amplitude of the residual 
correlation by 10%-20%. The PSF solution on each CCD 
on each exposure is iteratively refined by comparing the 
model PSF with the observed star and eliminating sig- 
nificant outliers. 

2.3.2. Step 2: PSF STACKING 

The PSF models for individual RESAMP images are 
inverse-variance weight-averaged to create the PSF on 
the mosaic image, where the weak-lensing signal is mea- 
sured. The image header of the RESAMP file contains the 
shift information (i.e., integer offsets). For each object, 
we need to loop over the list of the RESAMP files to stack 
PSFs. In order to determine whether or not an object is 
observed by a given RESAMP image and also to find the 
exact weight value used in co-adding, we utilize the corre- 
sponding "projected" weight map generated by Swarp. If 
the object is found to be within the weight map, we com- 
pute the PSF at the shifted location (remember that the 
RESAMP images are already rotated to properly align with 
the final stack) and applied the corresponding weight. 

Figure [5] and [6] illustrate that the above PSF stacking 
scheme closely reproduces the observed ellipticity pattern 
in F2. In Figure [5l the whiskers show the ellipticity dis- 
tribution of the stars directly measured from the 2° x 2° 
mosaic image. The PSF ellipticity change pattern mim- 
ics the 3x3 pointing pattern of the DLS observation. In 
addition, it is easy to see that across the exposure bound- 
aries (where the level of the shade changes) the PSF el- 
lipticity often exhibits a sudden change. The whiskers in 
Figure |6] display the ellipticity of the model PSFs evalu- 
ated (interpolation + stacking) at the location of these 
stars. The similarity in both the size and direction of the 
whiskers across the entire field is remarkable. 

Despite this seemingly nice agreement in the PSF on 
the stacked image, however, we find that this initial PSF 
model must be "tweaked" to remove the PSF-induced 
anisotropy to our satisfaction. This tweaking is carried 
out in two steps, for which we provide the details as 
follows. 

First, the model PSF tends to have systematically 



lower ellipticity by 5e = 0.001 ~ 0.003 with respect 
to the data PSF. This is because the procedure in the 
PSF sampling from noisy stars slightly circularizes the 
model PSF. Using this imperfect PSF model for our 
galaxy shape measurement leads to non-negligible under- 
correction. Hence, we compensate for this circularization 
by increasing the ellipticity (without altering the position 
angle) of the model PSF by 6e = 0.001 ~ 0.003. This 
"re-stretching" is implemented by shearing the PSF im- 
age in real space, and the applied shear is a constant 
(fixed for each DLS field) fraction of the PSF ellipticity. 

The exact amount of re-stretching for this first-level 
tweaking is determined using the following two diagnostic 
functions proposed by Rowe (2010): 

D^{r) = ((e<j - e„)*(ed - e„)) (r) (2) 
D2{r) = (e*(ed - e„0 + (e^ - e™)*ed) (r) (3) 

where and em are the ellipticity of the data and model 
PSFs, respectively in complex notation (see ij2.5.1|) . Con- 
sequently, Di and D2 show the residual autocorrela- 
tion and the data-residual cross-correlation, respectively. 
Rowe (2010) suggests that a combined use of these two 
functions provides an insight into systematics of the 
model. In Figure [7] we display Di{r) and D2{r) for F2. 
The left panel displays the result directly obtained from 
our PSF stacking whereas the middle panel shows the re- 
sult when this PSF model on the left panel is re-stretched 
to compensate for the PSF circularization. The improve- 
ment is more noticeable in D2{r). When comparing the 
amplitudes of Di and 1)2, one should remember that 
D2 (r) is in general more sensitive to the presence of sys- 
tematics than Di{r) in part because D2{r) is a sum of 
two data-residual ellipticity correlation functions (in or- 
der to cancel the imaginary part), and in part because 
the ellipticity of the PSF is higher than that of the resid- 
ual. For other possible reasons, we refer readers to Rowe 
(2010). 

The small residual correlation functions in the middle 
panel of Figure [7] suggests that the above re-stretched 
PSF model is an excellent description of the data. How- 
ever, we notice that the shapes of galaxies obtained with 
this re-stretched PSF (middle) tends to be still under- 
corrected. In other words, collectively speaking, galaxy 
shapes are still biased toward the initial anisotropy of 
the PSF. We suspect that this phenomenon is in part 
related to the so-called centroid bias mentioned by Bern- 
stein & Jarvis (2002) and Kaiser (2010). According to 
these authors, even a perfect PSF model will not remove 
the PSF bias completely because the centroid of the ob- 
ject is more uncertain along the elongation of the PSF. 
This bias does not go away even if we treat the centroid 
as free parameters. This is the reason that we need a 
second-level tweaking mentioned above. We address this 
issue by further stretching the model PSF so that the 
ellipticity increases by additional {5e) ~ 3 x 10"''. The 
amount of this additional stretching is also a fixed (for 
each DLS field) fraction of the PSF ellipticity, and the 
first-order value is determined mainly utilizing our im- 
age simulations, where galaxies are randomly oriented 
(i.e., no shear is present). We adjust the stretching fac- 
tor until the PSF-induced residual shear signal vanishes. 
Then, we refine this factor by making sure that the am- 
plitude of star-galaxy correlations ( ii3.2.1|) and B-mode 
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Figure 5. Observed PSF ellipticity in the stacked image for F2. The whiskers display the ellipticity distribution of the stars directly 
measured from the 2° X 2° mosaic image. The background shade represents the weight map (darker shade indicates lower value) derived 
from both exposure maps and photon statistics, and illustrates the complexity of the weight distribution. PSF discontinuities occur at 
exposure boundaries (i.e., at discontinuities in the weight map). We did not clean up outliers (e.g., cosmic-ray hit stars, binary stars, etc.), 
and they do not represent real PSFs. 



signals ( §3.2.2|) also decreases simultaneously. The right 
panel of Figure [7] shows the resulting Di (r) and D2 (r) 
diagnostic functions when this second-tweak is applied 
to the first-tweak PSF model shown in the middle panel. 
Note that this increases the deviation of D2{r) from zero 
at 8' < 6* < 70', ahhough this final PSF removes the 
PSF-induced anisotropy from galaxy images most satis- 
factorily among the three cases shown here. Figure [8] 
displays the star-galaxy correlation functions (see i i3.2.1l 
for the definitions) for the three cases shown in Figure [71 



It is obvious that the PSF model that we obtain from 
the second tweak gives the smallest amplitude for star- 
galaxy correlations, although the amplitude of the di- 
agnostic function (especially D2) of this PSF is not the 
smallest. 

One should not be misled into thinking that our PSF- 
tweaking removes any arbitrary B-mode signal. System- 
atics arising from non-centroid bias cannot be made dis- 
appear by simply increasing the ellipticity of every model 
PSF uniformly by a constant factor. In addition, the 
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Figure 6. PCA PSF reconstruction in the stacked image for F2. The whiskers show the eUipticity of the PSFs evaluated (interpolation 
+ stacking) at the location of the stars in Figure \E\ The agreement between observation and model in both the size and direction of the 
whiskers across the entire field is remarkable. 



above PSF-tweaking cannot arbitrarily get rid of intrin- 
sic alignment signals. 

2.4. Galaxy EUipticity Measurement 

There exist a number of algorithms for galaxy shape 
measurement in the context of weak lensing. Depend- 
ing on how one approaches the issue for the removal 
of PSF effect, we can classify the existing algorithms 
into moments-based methods and profile-fitting meth- 
ods. The former methods measure second-moments for 
both galaxies and PSFs and use them to estimate the 
pre-seeing eUipticity. This approach was pioneered by 



Kaiser, Squires, & Broadhurst (1995) and Fischer & 
Tyson (1997), and many variations exist. The later al- 
gorithms approximate the surface profile of galaxies with 
some analytic profiles. These analytic profiles are con- 
volved with PSF models before being fit to the images 
rather than fit to a deconvolved image. While the clas- 
sic, moments-based methods continue to be popular and 
updated, cosmic shear studies are relying more on the 
second, profile-fitting approach to overcome the poten- 
tial limitations (Kaiser 2000) of the moments-based ap- 
proach. 

Our shape measurement algorithm belongs to the sec- 
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Figure 7. Diagnostic ellipticity correlation functions for PSF modeling. -Di(r) and D2{r) are the residual autocorrelation and the data- 
residual cross-correlation, respectively (Rowe 2010). Any significant departure i^rom zero indicates that the model possesses non- negligible 
systematics. Here we show the case for F2. The left panel displays the result directly obtained from our raw PSF model whereas the 
middle and right panels shows the results obtained after the application of the first-level and second-level corrections, respectively. The 
first tweak is needed to improve the PSF ellipticity agreement between the model and data. This is done by increasing the ellipticity of the 
model PSF by 0.001 — 0.002. However, this first-tweak model does not remove the PSF-induced anisotropy in galaxy images completely 
due to the centroid bias. We have to further increase the ellipticity of the first-tweaked PSF by ~3 X 10"'* to remove this residual bias; 
note that this makes D2{r) deviates from zero at 8' < 9 < 70'. One should remember that D2{r) is in general more sensitive to systematics 
than Di{r) in part because the intrinsic ellipticity of the PSF is much larger than the residual PSF ellipticity (see Rowe 2010 for extensive 
discussion on the issue). 
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Figure 8. Star-galaxy correlation as diagnostics of PSF model. 
to th e PSF models shown in the left, middle, and right panels in i 
^3.1l for the definition) correlations, respectively. 

ond category. We fit a PSF- convolved elliptical Gaus- 
sian to a galaxy image. Of course, an elliptical Gaussian 
profile is not the best approximation of galaxy profiles. 
This sub-optimal fitting is termed "underfitting" (Bern- 
stein 2011) and has been shown to cause some bias in 
shear estimation. However, we find that this bias is only 
multiplicative and thus can be calibrated out with care- 
ful image simulations (discussed in t j2.5.2l) . Our experi- 
ments with Sersic profiles show that although this mul- 
tiplicative factor is reduced, the measurement uncertain- 
ties increase. This increase in ellipticity uncertainty is 
attributed to the following two facts. First, Sersic profile 
fitting takes into account more pixels farther from the 
object center, introducing larger noise. Second, Sersic 
profile fitting involves more free parameters to marginal- 
ize over. We want to include as many faint galaxies as 
possible for shear measurement as long as the net noise 
(quadratic sum of systematic and statistical noise) goes 
down, and we find that using Gaussian over other more 
sophisticated profiles increases the overall S /N of our cos- 
mic shear signal. 

Formally, a description of a galaxy image with an ellip- 
tical Gaussian requires the following seven free param- 
eters: normalization, semi-major and semi-minor axes, 
position angle, background level, and two parameters 



'e show the case for F2. The left, middle, and right panels correspond 
igure[7] Filled and open circles correspond to the "tt" and "xx" (see 

for the centroid. We fix the centroitfl and the back- 
ground using the SExtractor's xwin_image, ywin_image, 
and background so that the total number of free param- 
eters is only four, which further stabilizes the minimiza- 
tion and reduces the ellipticity uncertainty. The initial 
guesses for these four parameters are computed utilizing 
SExtractor measurements. 

For each object, square postage stamp images are ex- 
tracted from the final stack and rms map. We choose the 
size of this postage stamp image to be (8a -t- 20) pixels on 
a side, where a is the semi-major axis initially determined 
by SExtractor. In most cases, the image contains pixels 
belonging to other objects and we need to mask them 
out. This is implemented by replacing the rms values of 
these pixels with very large numbers, thus masking them 
out in further processing. The identification of these pix- 
els is based on the information in the segmentation map 
output by SExtractor. The shape measurement code is 
written in IDL, and the MPFItIJ module was employed 
as a minimizcr. MPFIT estimates parameter uncertain- 
ties from a Hessian matrix. We convert these errors to 
ellipticity uncertainties by error propagation. 

^ When we free the object centroid, the number of usable galaxies 
decrease by ~8% and the multiplicative shear calibration factor 

increases by 3 — 6%. 

^ available at |http:/ /www. physics. wisc.edu/ ~|craigm/idl/. 
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2.5. Shear Estimation 
2.5.1. Shear Estimator 

Gravitational lensing transforms the shape in the 
source plane to the image plane according to the fol- 
lowing matrix: 



(1 



1 - .91 -.92 
-.92 1 + .91 



(4) 



where k is the projected mass density in units of the 
critical lensing density and g is the reduced shear g = 
7/(1 — k). In the weak-lensing regime, k is small and 
thus the 7 — g assumption is often made. The (1 — k) 
factor affects the overall magnification, which is observ- 
able through the measurement of bias in object num- 
ber density or size distributions. The transformation 
matrix shears a circle into an ellipse with an elliptic- 

ity 9 = (.9i + 52)^^^ = (1 ~ + where r is 

the ratio of the semi-minor axis to the semi-major axis 
(i.e., b/a). The position angle of the ellipse is given by 

1/2 tan-i (52/51). 

Using complex notation g = gi + ig2, we can also ex- 
press the ellipticity transformation when an object has 
an initial ellipticity e = ei -I- ie2 as 



e = 



1 + g*e' 



(5) 



where the asterisk represents complex conjugation and e' 
is the measured ellipticity. If we assume that the distri- 
bution of e is isotropic, we can derive g from averaging 
over a population of galaxies using 



g = 



(6) 



where is a weight for each galaxy i. In the current 
paper, we use the following inverse variance as weight: 



1 



{6ei 



(7) 



where asN is a shape noise of the population per compo- 
nent (~0.25) and 6ei is the ellipticity measurement error 
per component. In equation [51 TZ is called the shear re- 
sponsivity, which is a calibration factor necessary to rec- 
oncile the difference between the average ellipticity and 
the shear. It is easy to show that 7^ « 1 if no mea- 
surement noise is present and galaxy morphology can be 
described by a simple elliptical isophote. However, be- 
cause neither is true in the real world, one must estimate 
TZ with care, and this is one of the most critical issues in 
future large lensing surveys since the result will not be 
limited by statistical uncertainties. 

Ideally, it is desirable to estimate TZ analytically from 
first principles and use image simulations only to ver- 
ify the accuracy. Bernstein & Jarvis (2002) provided an 
important contribution and their prescription has been 
used in quite a few studies (e.g., Jarvis et al. 2006, Hi- 
rata et al. 2004). Nevertheless, more rigorous efforts are 
still needed to quantify the applicability and limits of the 
Bernstein & Jarvis (2002) method. For the current DLS 
cosmic shear analysis, we find that the shear responsiv- 
ity derived with the Bernstein & Jarvis (2002) method 
agrees reasonably well with the value obtained from 



our weak-lensing image simulations for bright (R < 22) 
galaxies, but gradually underestimates the shear dilu- 
tion effect as the S/N of the objects decreases. Our DLS 
shear calibration hereafter is purely based on our image 
simulation studies, which are described in detail below. 

2.5.2. Image Simulation 

The translation of the measured ellipticity to the ap- 
plied shear is not straightforward. First, a response to 
a shear depends on galaxy populations. This is because 
the change in the second moments under a given shear 7 
depends not only on the second moments themselves, but 
also on the higher moments (Mandelbaum et al. 2012). 
This makes the effects of morphological features such as 
radial profiles, bulge-to-disk ratios, spiral arms, etc. non- 
negligible. As we model a galaxy light distribution with 
an elliptical Gaussian in the current study, we should 
understand how much the lack of details in the model 
biases the lensing signal. Second, ellipticity measure- 
ment is a noisy process. As most lensing signals come 
from faint galaxies, this measurement noise significantly 
dilutes the signal. Third, a nontrivial fraction of galaxies 
are affected by catastrophic shape measurement errors. 
The sources of these catastrophic shape errors include 
substructures of galaxies (e.g., HII regions), crowding, 
"bleed" trails, clipped objects, galactic cirrus, spurious 
detection around bright objects, etc. As the ellipticity 
measurement from these sources does not contain any 
lensing signal, the direction of the bias will always be to- 
ward underestimation. In the current paper, instead of 
quantifying the effect of each factor separately, we choose 
to derive a global value for shear responsivity TZ. Al- 
though it is worth investigating the effect of each factor in 
isolation, marginalizing over other parameters increases 
the required number of simulated image sets consider- 
ably, which is beyond the scope of the current study. 

We utilize a modified version of the Large Synoptic 
Survey Telescope (LSST) image simulator presented in 
Jee & Tyson (2011). The simulator samples galaxy im- 
ages from the Hubble Space Telescope (HST) / Ultra 
Deep Field (UDF; Beckwith et al. 2003) images and 
convolves them with the PSFs computed from the at- 
mospheric turbulence model and the telescope optics. 
The purpose of this modified image simulator is to cal- 
ibrate the conversion of ellipticity to shear. Given the 
same galaxy profile, the size and intrinsic ellipticity of 
the PSF are the most important factors affecting this 
calibration parameter. The main difference in the PSF 
between LSST and the two 4-m telescopes comes from 

1) different f-ratios (f/1.24 and f/2.7 for LSST and May- 
all, respectively), 2) exposure time (15 s vs. 900 s), and 
3) atmospheric seeing (0'.'65 versus 0'.'85). We address 

2) and 3) by changing the atmospheric parameters (e.g.. 
Fried parameter and outer scale) in such a way that the 
resulting seeing distribution is close to the observation. 
We cannot address 1) directly without replacing the cur- 
rent LSST optical design model with the most up-to-date 
Mayall/Blanco telescope models. However, it is possible 
to approximate the effect by degrading the focus (and 
optical alignment) so that when the diffraction limited 
PSF is convolved with the atmospheric PSF, it matches 
the DLS pattern. Without this adjustment, the delivered 
DLS PSF is severely circularized by atmosphere (longer 
exposure and large atmospheric PSF). After this modifi- 
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Figure 9. Shear recovery test. The calibration factor (inverse 
slope of the solid lines) depends on the S/N of the source galaxies. 
The high S/N galaxies (median magnitude ~ 23) require a calibra- 
tion factor of ~1.05 whereas the calibration factor of the low S/N 
galaxies (median magnitude ~ 26) is as high as ~ 1.28. 

cation, we obtain a distribution of PSF ellipticity ranging 
from 2% to 7%, matching the DLS data. Another im- 
portant question might be whether or not the resuhing 
spatial variation within a single DLS CCD is realistic. If 
our simulated PSF lacks a small scale variation compared 
to that of the data, the PSF model in the simulation may 
be easier to describe than in real situations. We find that 
the residual PSF correlation for both simulation and data 
shows a similar residual amplitude at all relevant scales, 
which suggests that the spatial variation of the PSF on 
both simulated images and DLS data possess a similar 
level of complexity. 

Most lensing signals come from z ^ 1 galaxies, the me- 
dian of the n{z) counts, which are dominated by the faint 
blue galaxy (FBG) population. Hence, it is important to 
test shear measurement from these UDF galaxies rather 
than from synthetic galaxies with analytic profiles. We 
randomize both the orientation and the position of the 
HST galaxies so that the net shear vanishes. We refer to 
these images as zero shear sky (ZSS). In a strict sense, the 
ZSS images are already convolved by the PSF of the Ad- 
vanced Camera for Surveys, and thus are not the "true" 
sky images in the absence of the instrument seeing. How- 
ever, because the size of the ACS PSF is a factor of eight 
smaller, the effect in the creation of the final DLS im- 
ages is limited to a scale far smaller than the DLS pixel 
0.257"). We apply a gravitational shear using bi-cubic 
interpolation to the ZSS images. As we have not down- 
sampled the UDF images yet, any interpolation artifacts 
and their propagation are expected to be insignificant on 
the final image. Then, we convolve these sheared sky 
(SS) images with spatially varying DLS PSFs. Readers 
are referred to Jee & Tyson (2011) for the details of the 
algorithm involved in this step. Finally, we down-sample 
the convolved sheared (CS) images and add noise to the 
result in order to match the pixel scale and depth of our 
DLS images. 

Our goal is to determine the relation between input 
shears and weighted sum of the ellipticities as a function 
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Figure 10. ij-band magnitude versus shear calibration factor. 
The data points ( diamond ) are measured from image simulations. 
The black solid line represents our parameterization of the shear 
calibration factor. The mean shear calibration factor derived by 
multiplying this curve and the weighted magnitude distribution of 
source galaxies (red) is ~1.08. 

of magnitude. Figure |9] shows the results for two of these 
simulated populations where the mean apparent magni- 
tudes are approximately 23 (open) and 26 (filled). We 
omit the results from intermediate magnitude objects to 
avoid clutter. The slope of the line is the shear responsiv- 
ity in equation [51 For the bright population, we obtain 
TZ ^ 0.95. This is similar to the value TZ ~ 0.93 that 
we obtain using eqn. 5.33 of Bernstein & Jarvis (2002). 
However, for the faint population the shear responsivity 
TZ ~ 0.78 determined from our image simulation is lower 
than the analytic estimate TZ ~ 0.89. Figure [10] summa- 
rizes the results when we combine the results of this shear 
recovery test for four magnitude bins. We parameterize 
the dependence of the multiplicative factor m-y — 1/7?. on 
the r— band magnitude (mj^) with the following form: 

= 6 X 10-'*(mi? - 20)3-2^ + 1.036, (8) 

where is SExtractor's MAG_AUTO. In deriving an aver- 
age multiplicative factor {m^), we need to consider both 
the magnitude and weight distributions of the source 
population; the source selection criteria are discussed 
in detail in i j2.7l The magnitude distribution of our 
source galaxies peaks at mn ~ 24 and then precipi- 
tously decreases (virtually no galaxies beyond 26). 
In addition, smaller weights are given to ellipticities of 
faint galaxies (eqn. [5] and [7]). The red histogram dis- 
plays this weighted magnitude distribution. The re- 
sulting mean multiplicative factor is estimated to be 
(m^) = 1.08 ±0.01. 

As this multiplicative factor is determined from the 
galaxies in the UDF, it is possible that the above shear 
calibration may need to be refined further for the DLS 
data. Indeed, our experience suggests that galaxy mor- 
phology and size distributions are non-negligible factors 
in shear calibration. Currently, the UDF images are the 
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only available space-based images that contain virtually 
noiseless galaxy images down to the limiting magnitude 
of the DLS. Therefore, to assess the effect of the sam- 
ple variance, we perform another suite of lensing image 
simulations, but this time with an analytic description 
of galaxy profiles. We utilize two publicly available soft- 
ware packages Stuff and SkyMaker, which create astro- 
nomical catalogs and images, respectively 0- The default 
generation of the ellipticity distribution in the Stuff cat- 
alog is rather unrealistic, and thus we modify the output 
in such a way that the ellipticity distribution per compo- 
nent matches the one in the UDF data. In addition, we 
force bulge and disk to align with an identical axis ratio. 
Gravitational lensing shear is applied at the catalog level 
by altering the object ellipticity. This allows us to mini- 
mize the dilution of the lensing signal from the interpola- 
tion noise. We first create space-based images and then 
convolve the results with the DLS-like PSF. The rest of 
the simulation follows the steps in our UDF-based anal- 
ysis. From this second set of simulations, we determine 
the mean multiplicative factor to be (m^) = 1.05 ± 0.01. 
As we observe that galaxies with analytic profiles tend to 
require smaller calibration factors, the ~3% decrease in 
(m^) in the latter experiment is consistent with our ex- 
pectation. Because it is unlikely that most of the galax- 
ies in the DLS can be decomposed into bulge and disk 
as are done here, we can assume that the difference in 
the galaxy population between the first and second sets 
of simulations represents an extreme case0. Therefore, 
we adopt the difference in (m-y) as the maximum devia- 
tion due to the sample variance, and we marginalize over 
1.05 < m~f < 1.11 with a flat prior in our cosmological 
parameter estimation. 

We did not participate directly in previous community 
shear calibration efforts such as STEP and GREAT, al- 
though it is worth mentioning here that our elliptical 
Gaussian fitting method is similar to the BJQ2 method 
where the authors propose to shear a galaxy image iter- 
atively until it matches a circular Gaussian. These shear 
calibration programs provided important contributions 
by raising the public awareness regarding the key issues 
of future weak-lensing surveys. However, our indepen- 
dent simulations include the following real- world features 
that both STEP and GREAT have not fully addressed 

yet.^ 

First, our training set data include a spatially varying 
PSF and its estimation through noisy stars. No spatially 
varying PSF was addressed in STEP and GREAT08 (Bri- 
dle et al. 2009). GREATIO (Kitching et al. 2012) ad- 
dressed the issue but only in a limited way. In the Galaxy 
Challenge where the participants are asked to measure 
pre-seeing ellipticity, the spatially varying PSF was given 
as a known function. In the separate challenge called 
the Star Challenge, the PSF at star positions was also 
provided as a known function. The only challenge in the 
latter is to interpolate/extrapolate the PSF to the galaxy 
location using the known PSFs at star positions. In real- 
world weak-lensing analysis the PSF must be estimated 
from a finite number of noisy stars, and this imperfect 

* http:/ /www. astromatic.net/software 

® When we relax the bulge-disk-aUgnment constraint, the mean 
multiplicative factor becomes {m-y) ~ 1.06, moving closer to the 
UDF case (m-y) ~ 1.08. 



model PSF is then applied to measure pre-seeing galaxy 
shapes. 

Second, we address the effect of galaxy morphology 
on shear measurement by using both real galaxy images 
and analytic profiles. As mentioned above, our simu- 
lation finds that at extreme cases the difference in the 
multiplicative factor is ~3% (i.e., 1.08 vs. 1.05). Since 
the DLS is deep, it is important to include many faint 
(> 24 ABmag) galaxies, whose high S/N proxy images 
are only available in the UDF data. 

Third, we simulate the effect of object blending. Both 
GREAT and STEP have assumed that a galaxy is iso- 
lated from the rest of the objects. However, when a sur- 
vey goes deep as the DLS, a significant fraction of the 
objects overlap with one another. Obviously, the de- 
tails of how one treats this blending of objects in source 
detection and shape measurement affect shear calibra- 
tion. Since we use realistic sky images (i.e., overlapping 
sources) and the same source detection/measurement al- 
gorithms for both DLS and training-set data, our shear 
calibration bias from blended objects is not an issue. 

2.6. Photometric Redshift Estimation 

Since the strength of gravitational lensing signal de- 
pends on the distance ratios between the observer, lens, 
and source, one's ability to characterize the amount of 
bias in photometric redshift estimation is as critical as 
the ability to control shear systematics for any preci- 
sion cosmic shear analysis. Because the lensing kernel 
is broad, individual galaxy photo- z estimation errors are 
of less concern than a skewed probability distribution. 
On the other hand, substantial effort should be made 
to address catastrophic errors, which can bias knowledge 
of the overall redshift distribution. The photometric red- 
shift catastrophic errors often arise when there are multi- 
ple peaks in the estimated probability distribution p{z). 
In many cases, these catastrophic errors are inevitable 
because of the inherent degeneracy between galaxy colors 
and redshift. In order to properly interpret the cosmic 
shear signal amplitude, it is important to understand the 
direction of the bias and quantify the fraction of catas- 
trophic outliers. 

We thus stack probability distributions of individual 
galaxies (instead of single-point, best- fit values) to re- 
construct the final redshift histogram of our source pop- 
ulation using the following equation: 

P(z) = A^P,(z)m. (9) 

i 

where pi (z) is the redshift probability distribution of an 
individual galaxy, fii is the weight used for our shear esti- 
mation (see §2.5.1|) . and A is the normalization constant. 
As argued by Wittman (2009), stacking p{z) provides a 
way to fairly represent the population with multimodal 
distribution in their p(z). In addition, even for galax- 
ies with unimodal redshift distribution, their p{zys are 
asymmetric in many cases because of the nonlinear map- 
ping of color space into redshift. Of course, one can con- 
sider discarding source galaxies with a multi-modal p{z) 
distribution or very broad probability peak. Although 
this may reduce the systematics effectively, the method 
decreases the number of galaxies containing gravitational 
lensing signals, resulting in a significant reduction in the 
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statistical power. 

We use the publicly available Bayesian Photometric 
Redshift code (BPZ; Benitez 2000). The six CWW+SB 
SED templates (Coleman et al. 1980; Kinney et al. 1996) 
enclosed with the BPZ code are tweaked so that we im- 
prove the agreement between best-fit values and known 
spectroscopic redshifts. We utilize the Smithsonian HEc- 
tospec Lensing Survey (SHELS; Geller et al. 2005) spec- 
troscopic redshift data (complete down to ^ 20.7) 
for this template tweaking. 

The advantage of the BPZ code is the use of magni- 
tude priors to partially break the color-redshift degener- 
acy We use the data from the VIMOS-VLT Deep Survey 
(VVDS; Le Fevre et al. 2005) to obtain magnitude- and 
type-dependent priors for our DLS photometric redshift 
estimation: 



p{z\m, T) = p{T\m) x p{z\T, m), 



(10) 



where the type dependence p(T\m) is parameterized for 
three types (E, Sp, and Im/SB) as 

p{T\m) = ft exp[-kt{m- 20)] (11) 

and the redshift dependence p(z|T, to) is parameterized 



p(z|m, T) — exp 



(12) 



In equation 1111 /( and kf are the type- dependent con- 
stants. Zm in equation 1121 is the median redshift for the 
magnitude to. Note that the faint tails of the priors are 
constrained strongly by the above functional forms, not 
by the small number of faint galaxies in the VVDS sam- 
ple. The comparison of our priors with those from the 
Hubble Deep Field (HDF) shows that the two sets of 
priors are very similar to each other. If we switched to 
HDF priors, this would shift the mean redshift of our 
DLS source galaxies by ~3%. 

As a spectroscopic test for fainter galaxies, we use an 
independent survey. The comparison between our DLS 
photo-z and the PRIsm MUhi-object Survey (PRIMUS; 
Coil et al. 2011) spec-z is shown in Figure [TT] for one 
field. The SHELS galaxies are ~2 mag brighter than the 
PRIMUS galaxies. Hence, this PRIMUS spec-z vs. our 
DLS photo-z comparison is a fair method to evaluate the 
performance of the above template tweaking carried out 
with the SHELS data. We compare 10,231 objects in 
F5, whose spectroscopic redshifts were measured by the 
PRIMUS project and quality flags (ZCDNF) are greater 
than 2. Although some systematic errors are indicated 
at the low {zspec ^ 0.15) redshift end, the overaU agree- 
ment is excellent between Zgpec and Zphot- The PRIMUS 
data is approximately complete down to to^ ^ 23 and 
about ~1, 100 redshifts are available for the objects in 
the 23 < rrir < 24 range. Hence, it is reasonable to 
assume that the DLS photo-z's are of good quality for 
the niR < 24 galaxies, which accounts for ~57% of our 
source population. We estimate that the aforementioned 
systematics at Zspec ^ 0.15 would affect the amplitude 
of the predicted shear signal by less than ^ 2% because 
the lensing kernel for our source population is broad. 
Beyond mn ~ 24, the fraction of catastrophic errors is 
expected to increase more rapidly with magnitude, and 
thus we argue that stacking p{z) should give a more re- 
alistic representation of the redshift distribution of the 



source population in this regime, too. We display the 
stacked p{z) of the entire source population in Figure fT2l 
Also plotted is the n{z) histogram computed from single- 
point photometric redshifts. The single-point photo-z 
histogram virtually truncates at z ~ 2, and this fea- 
ture is rather unphysical. Typically, cosmic shear studies 
parameterize the redshift distribution with an analytic 
form (e.g., p{z) cx (z/zO)" exp[— (z/zq)^]). Applying the 
method would smooth the distribution and recover the 
high-redshift tail as shown by the p{z) curve. 

The role of the magnitude prior becomes more im- 
portant progressively with magnitude as the SED con- 
straint weakens. At the faint end, we expect that the 
redshift probability p{z) of a nontrivial fraction of galax- 
ies may default to the prior. If we had taken a prior 
from small-field results such as the HUDF studies (Coe 
et al. 2006), the sample variance would have been a 
dominant source of bias for the rriR > 24 population. 
However, because the VVDS prior used here is obtained 
from a relatively large survey (~2,000 times larger than 
the HDFN in area), the impact of the sample variance 
on our photometric redshift estimation in DLS should be 
small. 

Therefore, it is fair to argue that the accuracy of the 
current p(z) for the to/j > 24 sources is limited by the 
systematics in the VVDS prior itself. Unfortunately, 
there is currently no solid method available for testing 
the systematics of the prior. Possible sources of system- 
atics include the effects of the galaxy SED evolution, 
inclination-dependent reddening, non-Gaussian photo- 
metric redshift errors, etc. For our cosmological parame- 
ter estimation, we assume a ~3% systematic error in the 
source mean redshift and marginalize over this interval; 
we implemented this by running a Monte Carlo Markov 
Chain (MCMC) while compressing/stretching the p{z) 
curve horizontally by a random factor drawn from the 
[0.97,1.03] interval for each chain. This 3% interval is 
the difference of the mean redshift when we compare the 
results from the VVDS priors and the HDF priors as 
mentioned above. In fact, the SED constraint is strong 
for a significant fraction of our source galaxies. There- 
fore, the above ~3% systematic error is a conservative 
value. 

2.7. Source Galaxy Selection 

Our source selection is based on various parameters. 
Table 2 summarizes the names of these parameters and 
the values used as threshold. The lower magnitude cut 
is needed to prevent an accidental inclusion of very low 
redshift galaxies and saturated objects. The upper limits 
for the magnitude and ellipticity error disallow spurious 
detections to be included. We set the lower bound of the 
photometric redshift to be Zphot > 0.3 because there is 
some indication that our systematics in photo-z estima- 
tion might be slightly biased at Zp^ot ^ 0.2, and because 
throwing away these low redshift galaxies minimizes a 
potential contamination of our cosmic shear from intrin- 
sic alignments. The lower limit in galaxy size is required 
to avoid stars and also very small galaxies whose shape 
measurements are noisy and biased because of the pixel- 
lation effect. In many cases, despite their misleadingly 
small ellipticity measurement errors, the ellipticities of 
these small galaxies show distinctive patterns as shown in 
Figure fT3l We threshold the semi-minor axis b measured 
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PRIMUS Spec-z Redshift 



Figure 11. Comparison of spectroscopic redshift with photometric redshift in DLS. In this plot, we display 10,231 objects in F5, whose 
spectroscopic redshifts were measured by the PRIMUS project and quality flags (ZCQNF) are greater than 2. The plot on the left panel 
compares the PRIMUS spec-^'s with DLS photo-z's directly whereas the plot on the right compares the distribution of the redshifts. The 
p(z) curve was obtained by stacking the photometric redshift probability distribution of individual galaxies. 



Table 2 

Parameters Used for Source Selection 



Name Condition 



_R-band magnitude 
ellipticity measurement error 
photometric redshift 
semi- minor axis 

elliptical gaussian fitting convergence 

source extraction flag 

masking 



21 < mfl < 27 

(5e < 0.25 

^phot > 0-3 

b > 0.4 pixels 

STATUS = 1 

FLAG ^ 8, 16, 32 

excluded if near or within the region 



from our elliptical gaussian fitting to be greater than -^0.3 
pixels in order to make the pixellation effect negligible. 
The MPFIT minimization engine reports the status of the 
convergence. Although the author comments that all 
values greater than zero in the keyword STATUS can rep- 
resent success, we select the objects only with STATUS 
= 1, which is the most conservative indicator of the con- 
vergence. The SExtractor software saves the history of 
the source extraction in the binary switch format in the 
FLAG parameter. We discard an object if any of the third 
(8), fourth (16), or fifth (32) bits is turned on. This helps 
us to exclude the objects on or close to image borders. 
We mask out the regions affected by bright stars (PSF 
wings and bleeding streaks). The resulting total number 
of source galaxies in our ~20 sq. degree area is ~1.2 mil- 
lion, which gives a mean number density of Usource ~ 17 
galaxies per sq. arcmin. Without the Zphot > 0.3 thresh- 
old, this number density would increase to risource 21 
galaxies per sq. arcmin. The field-to-field variation of 
Ubg is listed in Table 1. The magnitude distribution of 
the source population is shown in Figure 

2.8. Intrinsic Alignment and Luminous Red Galaxies 



A potentially important systematic in cosmic shear of 
astrophysical origin is intrinsic alignment of source galax- 
ies (e.g., Hirata & Seljak 2004). If a large-scale tidal grav- 
itational field significantly affects intrinsic alignments, 
our interpretation of cosmic shear measurements must 
include these intrinsic ellipticity correlations. Ellipticity 
correlations between galaxies subject to a common large 
scale gravitational field are often called the intrinsic- 
intrinsic (II) signal whereas the correlation between fore- 
ground galaxy ellipticity (tidal) and background galaxy 
shear (gravitational lensing) is called the gravitational- 
intrinsic (GI) signal. 

In tomographic studies, the auto-correlation function 
within a narrow redshift shell may be severely influ- 
enced by the II signal. However, in the current non- 
tomographic study, the II signal is expected to be negli- 
gible because galaxy pairs within a close redshift interval 
are substantially outnumbered by those with very differ- 
ent redshifts. In addition, because the II signal increases 
for decreasing angular scale, we can mitigate this po- 
tentially small contribution further by discarding close 
pairs. 

However, the GI signal may be non-negligible even 
in non-tomographic studies because the fraction of 
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redshift 

Figure 12. Redshift distribution of source galaxies. The his- 
togram constructed from best-fit photo-z values is compared with 
the p{z) curve obtained by stacking individual probability distribu- 
tions. We discard the sources with Zphot < 0.3 because our system- 
atics in photo-z estimation might be slightly biased at Zphot J; 0-2, 
and because throwing away these low redshift galaxies removes a 
large fraction of luminous red galaxies ( ^2.71 . 

foreground-background pairs is certainly overwhelming. 
Mandelbaum et al. (2006) report a significant detec- 
tion of the GI signal from the Sloan Digital Sky Sur- 
vey (SDSS). They conclude that luminous red galaxies 
(LRGs) are the main sources of the intrinsic alignment 
and project that cosmic shear surveys at z ^ 1 may un- 
derestimate the linear amplitude of fluctuations as much 
as 20%. 

We address removing LRGs from our source catalog 
by identifying the population with our photometric red- 
shift catalog. We classify galaxies as LRGs whose SED 
template is consistent with the elliptical galaxy template 
and absolute magnitude is brighter than Mn < —20. We 
limit the maximum redshift to be Zphot < 0.9 because 
LRGs at higher redshift should provide negligible contri- 
bution to the GI signal. From this procedure, we discard 
about 5% (~ 61, 200 objects) of our original source galax- 
ies. 

A comparison of the cosmic shear signals with and 
without the LRGs shows that the shift in signal ampli- 
tude is much smaller than the statistical errors and thus 
does not affect our cosmological parameter estimation. 
This non-detection of the GI effect is rather unexpected. 
However, given the depth of the DLS and thus the fact 
that most lensing signals come from high redshift, it is 
possible that the GI contribution may become relatively 
insignificant. In addition, because we do not use galaxies 
at Zphot < 0.3, we suspect that the GI effect is already 
suppressed even before we remove LRGs. 

No significant detection of the GI effect was claimed 
in some previous studies. Fu et al. (2008) report that 
their marginalized likelihood analysis with the CFHT 
data does not show any negative GI signal as predicted 
by the theory. Schrabback et al. (2010) state that their 
redshift scaling analysis is not affected by the presence of 
LRGs. The early CFHT study by Fu et al. (2008) may 
have suffered from non-negligible weak-lensing systemat- 



ics (Heymans et al. 2012, Kilbinger et al. 2009). Also, 
the redshift scaling test by Schrabback et al. (2010) is 
not very sensitive to the GI signal. Hence, it is premature 
to draw any firm conclusion from these results. For the 
DLS, we have yet to perform a full analysis of the intrin- 
sic alignment, including source galaxies at Zphot < 0.3. 
We thus drop these galaxies from our sample. 

3. COSMIC SHEAR MEASUREMENT 
3.1. Theoretical Background 

We carry out cosmic shear analysis with two-point 
shear-shear correlation functions. These second-order 
statistics and their derived quantities have been widely 
used, and the mathematical tools and the algorithms are 
in a mature stage. The computation of correlation func- 
tions is time-consuming if brute-force algorithms are em- 
ployed. Therefore, we developed a fast tree-code, which 
closely approximates the result from the exact brute- 
force method. Our tree-code is cross-checked against 
some publicly available codes EB- We note that although 
we used our fast tree-code for intermediate steps, the fi- 
nal results presented in this paper are obtained through 
our brute-force two-point correlation estimation code. 

The shear-shear correlations are evaluated as 

Ue) = ^''^"^"^^^''^^-^ (13) 

and 

X {6) = ^».j"^»^'jex,»ex,j-^ ^^^^ 

I^WiWj 

where the summation is carried out over every possible 
combination of i*'* and j*'' galaxies {i < j). The two 
subscripts t and x refer to the two projections of the 
ellipticity along the tangential and 45° angle with respect 
to the line connecting the galaxy pair, respectively. 

is the weight associated with the ellipticity of the i{jY^ 
galaxy. The angle between the galaxy pair is 6. 

The two following linear combinations of ^tt and x 
are useful to express the derived shear statistics that we 
discuss hereafter: 

e+-6i+exx (15) 

and 

?-=6t-exx. (16) 
Derived shear statistics commonly used include shear 
variance (^Ie/b^^ aperture mass {^M'e/bJ^ and correla- 
tion function ^e/b, where the subscripts E and B de- 
note the contributions from the so-called E- and B-mode 
signals, respectively. The E/B signals are analogous to 
electric versus magnetic fields in the electromagnetic the- 
ory, where the former is the gradient of a scalar field and 
the latter is the curl of a vector field. Because gravita- 
tional lensing only produces E-mode (curl free) signals, 
in principle the B-mode signal must be consistent with 
zero. As image distortions due to systematics (e.g., in- 
accurate PSF correction) are not always curl-free, the 
B-mode measurement is frequently used as a diagnostic 
of the residual systematic errors. 

e.g.. |http: //ww w2.iap.fr/users/kilbinge/athena7| 

|http: / / code.google.eom/p/mjarvis7] etc. 
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Figure 13. Effect of object size on ellipticity measurement. Extremely small objects show distinctive patterns regardless of their S/N in 
their ellipticity distribution mainly due to the pixellation, which disappears for relatively large objects. We threshold the semi-minor axis 
b measured from our elliptical Gaussian fitting to be greater than 0.4 pixels in order to make the pixellation effect negligible. The objects 
in both panels have a "reported" ellipticity measurement error less than 0.25. 

The difference between equations [IT] and [18] is the use 
of different filter functions (5'+._ and T^ _), which sam- 
ple different angular parts of ^+,-. The filter functions 
S+,- and T_|___ are defined in Schneider et al. (2002). 
The E- and I3-mode decomposition for the correlation 
function requires the definition of the following quan- 
tity C': 
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Figure 14. ii-band magnitude distribution of source galaxies. 
We select a total of ~1.2 million galaxies from the five DLS fields. 
When we convert galaxy ellipticities into shears, the measured el- 
lipticity is weighted with the inverse variance, which is discussed 
in the text. The red histogram shows the magnitude distribution 
when this weighting is considered. The black dotted histogram 
displays the raw magnitude distribution without this weighting. 
The mean shear calibration factor should be derived using the red 
histogram. 

The top-hat shear variance is given by 

, , , /i3 



,{e) : 



26^ lo 



The aperture mass statistic is 



(17) 



E/I 



(S), 



1 n 

2e^ Jo 



«'(«) = ?-(e) + 4 



r 



129- 



i: 



1?3 ' 



which is also the result obtained by filtering 
<^', we can define £,e/b follows: 



^E,B{e) = 



(19) 

Using 
(20) 



The above three lensing statistics assume that we can 
measure the C-f/- data on arbitrarily small and large 
scales. Therefore, when these statistics are evaluated us- 
ing real data with a cutoff at both ends, the resulting 
E/B-decomposition values deviate from the theoretical 
ones (e.g., Kilbinger et al. 2006). This limitation moti- 
vates some authors to suggest new E/B-mode statistics 
that do not suffer from the finite-interval problem (e.g., 
Schneider & Kilbinger 2007; Eifier et al. 2010; Fu & Kil- 
binger 2010; Schneider et al. 2010; Becker 2012). Among 
these, we consider the ring statistics (Schneider & Kil- 
binger 2007) using a scale-dependent integration limit rj 
suggested by Eifier et al. (2010). The parameter refers 
to the ratio of the smallest separation to the largest sep- 
aration 77 = Omin/Q- The ring statistics is evaluated in a 
similar fashion as above except that the integration limit 
is over a finite interval [77^, 0] as the following; 



1 r" M 

{RRE/B)(e)=- —[Z+(§,r,)(,+ (§)±Z^(-d,ri)£,^(S)], (21) 
^ .1 r,e " 



(18) 



where the filter function _ is defined in Schneider 
Kilbinger (2007). 
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Now, in order to compare the above statistics obtained 
from our DLS data with the prediction for a given cos- 
mology, we need a method to predict the signaL We be- 
gin with the foUowing cosmic convergence/shear power 
spectrum: 

where Hq is the Hubble parameter, fli^i is today's matter 
density, a{x) is the scale factor at a redshift correspond- 
ing to a comoving distance x, Jk is the comoving angular 
diameter distance, Ps is the linear power spectrum, and 
gk is the lensing efficiency factor: 

for the fc*'' redshift shell with a redshift distribution 
Pkix)- Note that only the redshift probability Pk{x) con- 
nects the survey to the shear power spectrum. Also keep 
in mind that in the current non-tomographic study only 
a single broad redshift shell is used. 

With this convergence power spectrum F^', we can 
obtain predicted cosmic shear statistics for any of the 
above. For example, the basic correlation function ^+,- 
is the convolution of the shear power spectrum with a 
Bessel function Jo,4 (Crittenden et al. 2002). That is, 

1 f°° 

^+-i0) = Tr dUJo,4mP!:\e)- (24) 

In principle, the aperture mass statistics {j^l^/g^ and 

{^"e/b^ can be evaluated by substituting these 

functions into Equations [IT] and [181 respectively. How- 
ever, it is more convenient to obtain them by directly con- 
volving the shear power spectrum P^'' with correspond- 
ing kernels in a similar way to Equation 1241 

3.2. Weak-lensing Systematics in DLS 

One, if not the most, critical source of systematics in 
weak lensing is the inaccurate removal of the PSF effects 
from shear measurements. The problem arises either be- 
cause the PSF model is in error or because the removal 
procedure is suboptimal despite the correct model. In 
any case, the effects of this PSF correction error can 
be classified into two kinds: shear calibration and resid- 
ual anisotropy. The degree of the shear calibration and 
anisotropy systematics are often parameterized by the 
following multiplicative and additive factors c-y, re- 
spectively: 

j — m^j + c^, (25) 

where 7 and 7 are the weight-averaged ellipticity and true 
shear, respectively. The multiplicative factor can be 
also affected by the degree of complexity in galaxy shape 
modeling, but we do not distinguish the effect of PSF 
from that of the galaxy shape modeling in this study. 
Instead, we determine the global value of m-y from the 
image simulations described in i )2.5.2l 

3.2.1. Star-Galaxy Correlation 

The additive error is mostly caused by imperfect PSF- 
induced anisotropy removal. A useful diagnostic of the 



additive error is the following star-galaxy correlation 
(Bacon et al. 2003): 

^^y^^l^)' ^^^^ 

where e* is the ellipticity of uncorrected stars. The in- 
trinsic size of the PSF ellipticity depends on many con- 
ditions and varies widely. Equation makes sure that 
the amplitude of the star-galaxy correlation (e*7) is nor- 
malized by this intrinsic size of the PSF ellipticity to en- 
able a fair comparison with other observations. However, 
occasionally the intrinsic PSF correlation (e*e*) crosses 
zero where the transitions between positive and negative 
correlations occur. In these cases, ^sys can rise abruptly 
because of the small denominator if (6*7) is indeed uncor- 
related with (e*e*). We use a mean amplitude of (e*e*) 
in the 10' < 9 < 100' range to avoid this artifact. 

In Figure [TSl we display the star-galaxy correlation for 
all five DLS fields. Overall, the amplitude of the star- 
galaxy correlation is far smaller than the cosmic shear 
errors evaluated for the field. Note that at 6* ~ 30', 
the star-galaxy correlation amplitude in F5 is compara- 
ble to the statistical errors. However, the sample vari- 
ance error is still much larger even in this rare case. It 
is important to remember that the galaxy shapes used 
here are measured with the second-tweak PSF discussed 
in tj2.3.2l These small star-galaxy correlations strongly 
suggest that the residual systematics after our PSF cor- 
rection is not a concern in this DLS cosmic shear study. 

3.2.2. E- and B-mode Decomposition 

As mentioned in §3.1) the B-mode statistics provide 
useful insights into the residual systematics in gravita- 
tional lensing. However, the evaluation of equations 1171 
1181 and l20l formally requires the measurement of on 
large and/or small scales inaccessible to the DLS data. 
This is a well-known problem in cosmic shear studies. 
According to Kilbinger et al. (2006), the missing data on 
large scales make the ^e/b statistics valid only up to an 
offset over the entire range whereas those on small scales 

suppress the (^M^^^^ statistics on small scales. The ex- 
act amount of deviation with respect to the result from 
the ideal case depends on the cutoff angle, redshift distri- 
bution, and cosmology. Using the parameters specific to 
our DLS study, we examine these effects quantitatively. 

The missing data on small scales is relevaniFI to (^M'^^g^ 

and displayed in the left panel of Figure [HI Although 
our cutoff angle happens at 9min ~ 0.'15, the suppres- 
sion effect is visible out to 61 ~ 4'. This implies that 
our (observed) E/B-decomposition using this statistic is 
meaningful only at 6* > 4'. On the other hand, the £,eib 
statistics is affected by the missing data on large scales 
[9 > 100") and the effect is manifested as a constant off- 
set over the entire scale (right panel of Figure [T6| . This 
predicted offset in is -4.4 x 10"^ for the WMAP7 
cosmology, which is in excellent agreement with the ob- 
served value ~4.5 x 10~® when we carry out the E-/B- 
decomposition with the finite-field DLS data. Although 

Strictly speaking, the upper limit of the integration is 29, and 
this causes a slight deviation near the maximum angular separation 
because of the missing data between and 29. 
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Figure 15. Star-galaxy correlation. The diamond and triangle symbols represent the tangential ((£^7+)) and cross ((exTx )) components 
of star-galaxy pairs divided by auto-correlation functions of stellar ellipticity. Overall, the amplitude of the star-galaxy correlation is far 
smaller than the cosmic shear errors evaluated for the field. The sample variance is estimated using the Sato et al. (2011) method for 4 sq. 
deg. The statistical error is computed by counting the number of pairs in each bin. 
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Figure 16. Effects of missing data on the f g/g, \1e/b/ ^"^^ \^'^e/b/ statistics. Wc assume the WMAP7 cosmology and the source 
redshift distribution shown in Figure [T2I The left plot illustrates the effect on (^J^I^/B^ '^^'^ *° missing data on small scales {6 < 0fl5) 
whereas the right plot shows the effect of the missing data on large scales {6 > 100') on the ^ g/g statistics. The tophat shear variance 
{t^b/s) (middle) suffers from the missing data on both ends. The solid and dashed lines represent the E- and B-mode signals, respectively. 



the exact value of the offset depends on the assumed cos- 
mology, the variation within trg = 0.7 — 0.9 is compara- 
ble to the sample variance of the DLS. Therefore, in our 
presentation of the ^e/b statistics hereafter we show the 
results obtained by supplementing ^_ at ^ > 100' with 
theoretical predictions for the given redshift distribution. 
The addition of the synthetic data at ^? > 100' precisely 
cancels the ^iA x 10^^ offset. This is justified because 
the result is virtually cosmology- independent. However, 



when it comes to 



), since the predicted ^+^~ sig- 
nal at 9 < 0.'15 is sensitive to cosmology, we present the 
results without filling in the missing data on this small 

scale. Consequently, our DLS ^-^J;/^) values are sup- 
pressed at 6* < 4'. The tophat shear variance (^^^/g 
(middle panel of Figure [TBI) suffers from the missing data 
on both ends. For the evaluation of (j^/g), we only 



18 



Jee et al. 



supplement ^_|_._ at 6* > 100' with theoretical values. We 
note that the amount of the observed signal suppression 

on small scales for { M'i, , „ ) and ii%ir,i are consistent 



with the theoretical prediction. 

We present these cosmic shear statistics in Figure [T71 
where the DLS field number (F1-F5) runs from top to 
bottom. The displayed error bars represent only shear 
shot noise. The cosmic shear signal is clearly seen in 
all three ((7I;), (-^J;), ^_e) statistics whereas the cor- 
responding B-modes are all close to zero. The shear 
variance (\7|;)) measures the mean dispersion within an 
aperture and hence the data points are highly correlated. 
The signal shape from all five fields is similar and the 
amplitude is roughly proportional to the amount of large 
scale structure in each field. For example, (7I;) in F2 is 
nearly a factor of two higher at ^ 10' than the signal 
in the other fields, and this is consistent with the rather 
unusually large number of structures in F2 seen by both 
red sequence distribution and convergence map (Kubo et 
al. 2009). These structures are known to be at 2; 0.3 
and 0.6, and their angular scales are consistent with the 
--10' scale of excess shear correlation. Because we do not 
exclude the data in F2 in our cosmological parameter es- 
timation, it is possible that the results are slightly shifted 
toward high normalization, although the DLS field were 
randomly chosen. The signal from S^e statistics is similar 
to (7I;) except that the data points are less correlated. 
As expected, the field to field variation in signal is the 
largest in aperture mass (M|;), and the correlation be- 
tween the data points is the least among the three statis- 



tics shown here. Therefore, although the S/N of 

for each data point is lower than the other two statistics, 

the information content as a whole is in fact comparable. 

In order to maximize the constraining power of our 
DLS data on cosmology, we need to combine the mea- 
surements ^+ and ^_ from all five fields with a careful 
choice of weighting scheme. The errors are dominated 
by statistical errors on small scales whereas the sample 
variance is more important on large scales. For our DLS 
data, the transition of this error dominance occurs at 
~ 2' (see Figure [T5|) . We weight the measurements in 
each field with the following inverse variance: 



1 



w 



(27) 



'sv 



shot 



where asv is the sample variance, ashot is the shot noise, 
and fJsys is the residual systematic error estimate. The 
sample variance usv is a function of both cosmology and 
angle, and we discuss the issue in detail in i j4.2l The 



shot noise ashot is given by (T|^iV~°-^, where N is the 
number of all galaxy pairs used per angular bin. The 
residual systematic error is estimated from the analysis 
of B-mode signals and star-galaxy correlation. Although 
our residual systematics are negligible compared to sta- 
tistical errors and the sample variance, the star-galaxy 
correlation functions (see Figure [T5|) indicate that the 
star-galaxy correlation might be non-zero at small angu- 
lar scales. We model the trend with: 



asys{e) = iy~io^''{e/i')- 



(28) 



In Figure [18] we show the combined ^+ and ^_ cor- 
relation functions from all five fields. We assume the 



WMAP7 cosmology for the estimation of the sample 
variance asv here, which however is allowed to vary 
in our cosmological parameter estimation (Q. Fig- 
ures [T21 Uni and [5T] display the three derived statistics 

{^e/b)^ (^^b/s)' ^E/B, respectively. As men- 
tioned above, the missing data on a small scale {6 < 0.'15) 
suppresses and (m]^^^ on small scales, consis- 

tent with the prediction (Figure [H]). In addition, no 
constant offset is observed in (je/b') ^^'^ ^e/b because 
we supplement the DLS data with the synthetic data at 
9 > 100'. Finally, we display the ring statistic (^RRe/b) 
in Figure [211 The plot confirms that the B-mode sig- 
nal is also consistent with zero in this statistic that does 
not require any synthetic data. Together with the star- 
galaxy correlation presented in i )3 2.11 the absence of 
the B-mode signals here supports our success controlling 
wcak-lcnsing systematics. 

Note that in Figures [TS j [TO l [^ m\ and [H we scale 
the cosmic shear signal using the shear calibration result 
ft j2.5.2[) . The various lines represent the predicted signal 
for different as values while we assume the WMAP7 cos- 
mology for the rest of the parameters. The comparison 
of these predicted values with the data points indicates 
that the DLS data may favor a slightly higher normal- 
ization than the WMAP7 prediction (i.e., erg — 0.81). In 
particular, this tendency is conspicuous for the ^+ data 
points in the 3' < 6' < 20' range (Figure [T8)l . We do not 
find any evidence that this is caused by residual system- 
atics. Rather, we believe that the rich structures in F2 
(see Figure [T7|) is responsible for this signal excess. 

In deriving cosmological parameters from the DLS 
data, we use only the two ^+ and ^_ correlation func- 
tions, which are directly measured from our shape cata- 
log. Although it is possible to obtain our parameter con- 
straints based on pure E-mode statistics such as (7I;), 
(Af^), £_E, and (RRe), most theoretical studies so far 
nave focused on the correct evaluation of the covariance 
for ^-|_ and ^_ rather than for those derived statistics. 
In principle, a constant shear field can lead to a result, 
where ^+ = const and ^_ = while those E/B-niode 
statistics are unaffected (Schneider et al. 2010). How- 
ever, we are unable to find any evidence indicating that 
a constant shear field from residual systematics might 
be present in our DLS data (e.g., the star-galaxy ellip- 
ticity correlation should reveal this type of systematics). 
Hence, we believe that as long as the covariance matrix is 
constructed from robust error propagation, our param- 
eter constraints should yield virtually indistinguishable 
results regardless of the choice in lensing statistics. 

4. COSMOLOGY WITH DLS COSMIC SHEAR 

4.1. Shear Power Spectrum Model 

Cosmological parameter estimation is performed by 
comparing observed signals with those predicted by full 
simulation of the survey in different cosmologies. As dis- 
cussed in H3.ll the key component in model predictions 
is the shear power spectrum i-l^, which is obtained by 
integrating the matter power spectrum Pg along the line 
of sight weighted by lensing efficiency. Using a nonlinear 
matter power spectrum Ps is critical, and also it is well 
known that the details in the method for the computation 
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Figure 17. Cosmic shear statistics for each DLS field. The plots in the first, second, and third columns show shear variance, aperture 
mass, and correlation functions, respectively. Filled and open circles represent E- and B-mode signals, respectively. The error bars show 
only shot noise. No shear calibration has been applied yet. 
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Figure 18. correlation functions from the combined (Fl- 

F5) cosmic shear data. Filled and open circles represent 5+ and 
5_ signals, respectively. Both 5+ and error bars include the 
shot noise and the sample variance. Dashed (dotted) lines repre- 
sent predicted ^-i- (§— ) signal for different erg values while we as- 
sume the WMAP7 cosmology for the rest of the parameters. Shear 
calibration has been applied. 



Figure 20. Aperture mass variance from the combined (F1-F5) 
cosmic shear data. Filled and open circles represent E- and B- 
mode signals, respectively. The error bars show only shot noise. 
Dashed lines represent predicted signal for different erg values while 
we assume the WMAP7 cosmology for the rest of the parameters. 
Shear calibration has been applied. 
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Figure 19. Top-hat shear variance from the combined (F1-F5) 
cosmic shear data. Filled and open circles represent E- and 13-mode 
signals, respectively. The error bars show only galaxy shot noise. 
Dashed lines represent predicted signal for different erg values while 
we assume the WMAP7 cosmology for the rest of the parameters. 
Shear calibration has been applied. 

of the nonlinear power spectrum produce non-negligible 
differences in the parameter estimation. For example, 
the use of Peacock and Dodds (1996) shifts the value of 
(jg high by ~5% with respect to the case when the Smith 
et al. (2003) method is used. In this paper, we use the 
modified transfer function of Eisenstein & Hu (1998) that 
includes Baryonic Acoustic Oscillations (BAO) features 
and the Smith et al. (2003) "HaloFit" nonlinear power 
spectrum. 



Figure 21. ^e/b statistics from the combined (F1-F5) cosmic 
shear data. Filled and open circles represent E- and B-mode sig- 
nals, respectively. The error bars include both the shot noise and 
the sample variance. Dashed lines represent predicted signal for 
different trg values while we assume the WMAP7 cosmology for 
the rest of the parameters. Shear calibration has been applied. 

Note that according to Hilbert et al. (2009) and Hcit- 
mann et al. (2010), HaloFit (Smith et al. 2003) undcr- 
predicts the matter power on small scales, and for COS- 
MOS cosmic shear data, using HaloFit prediction leads 
to an overestimation of (ig by about 5% (Schrabback et 
al. 2009). Another unresolved issue in our theoretical 
cosmic shear signal modeling is the effect of baryons on 
the matter power spectrum, which affects the power spec- 
trum on small scales {I > 1000) (e.g., Jing et al. 2006, 
Rudd et al. 2008). 
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Figure 22. (^RR^^g') statistics from the combined (F1-F5) cos- 
mic shear data. Filled and open circles represent E- and B-mode 
signals, respectively. The error bars here include only shot noise. 
Dashed lines represent predicted signal for different erg values while 
we assume the WMAP7 cosmology for the rest of the parameters. 
Shear calibration has been applied. 

The shear power spectrum obtained in this way as- 
sumes that the reduced shear g is equal to the shear 7 
fi j2.5.ip . which leads to underestimation of the amplitude 
on small scales. Kilbinger (2010) provides fitting formu- 
lae, which prescribe the amount of correction necessary 
to improve the accuracy of the shear power spectrum. 
As the fitting formulae are valid for the cosmological pa- 
rameter set near the fiducial fiat ACDM WMAP7-like 
values, we cannot apply the correction to our model for 
the full range of parameters. However, the correction 
is very small (e.g., 1% at = 0f5), and we find that 
this correction has virtually no effect on our cosmological 
parameter estimation on the scale of the noise. 

We fix the baryon density and spectral index to be 
rib = 0.046 and Ug — 0.96, respectively (Komatsu et al. 
2011). The Hubble parameter h is marginalized over the 
0.64 < h < 0.80 range with a flat prior, consistent with 
the Hubble Space Telescope Key Project (Freedman et al. 
2001). We only consider flat (fijv/ -I- ^^A = 1) universes. 

4.2. Data Vector and Covariance Estimation 
Our two-point correlation data vector ^ consists of two 
parts: ^_)_ and The likelihood function is given by 



oxp -0.5 U-€ c-Vc-c 



(29) 



(27r)"/2|C|i/2 

where ^ is the model prediction for a given set of cosmo- 
logical parameters, C is the covariance matrix, and n is 
the number of elements in the data vector 4. The deter- 
minant |C| in the normalization should not be ignored 
in studies such as this one, where we take into account 
the cosmology dependence. 
The structure of C is 

^^(cfc::,)' ™ 

where the sub-matrices and C__ are the covariance 




Figure 23. Covariance m atrix of data vector. We represent the 
matrix element of Equation 1301 using the linear scale color scheme 
shown on the top. Only the contribution from the sample vari- 
ance is displayed for the fiducial WMAP7 cosmology. The angular 
bins are logarithmically spaced with i = I and i = n being ~0f3 
and ~100', respectively. The sample variance is dominated by the 
nonlinear part of ^!|_ on small scales [6 < 10'). Also, note that 

cov[^^,^^] and cov[^^,5^] are smaller than cov[5^,5^] by an or- 
der of magnitude. 

matrices of ^_|_ and ^_ , respectively, and the off-diagonal 

sub-matrix C-| is the covariance between and ^_ . 

Because we use the cosmic shear signals in the range 
0.3' < 6 < 100' comprised of 30 data bins, the dimension 
of C is 60 X 60 in our study. 

The covariance C is decomposed as C = C„ + Cs + 
Cb + Ce, where C„ is the statistical noise, Cb is the 
residual systematic error, Cs is the sample variance, and 
Ce is a cross-term between shape noise and shear corre- 
lations (Joachimi et al. 2008). C„ is directly measured 
from our DLS shape catalog {cr^j^/N). We assume that 
Cs is a diagonal matrix whose elements are given by 
Equation 1281 The contribution from might become 
important in some cases. However, for our DLS study 
we find that including these cross terms only negligibly 
affects our cosmological parameter constraints. This can 

be understood because Ce mostly contributes to the C_| 

covariance between large-scale and small-scale ^+ val- 
ues. Since the constraining power of is insignificant in 
our DLS case, our cosmological contours virtually remain 
the same, regardless of the presence of these terms. 

The Gaussian components of can be easily derived 
from the shear power spectrum. For example, the Gaus- 
sian covariance of £,+ {9) and ^+(6*') is given by (Joachimi 
et al. 2008) 

CovK+(e),c+(e')l = f i di j„(w)Mw')p^{if, (31) 

where A is the sky area of the survey. Because the avail- 
able Fourier modes are limited by the area, equation |311 
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in fact leads to overestimation. According to Sato et 
al (2011), the discrepancy is roughly a factor of two for 
A ^ 25 deg . Sato et al. (2011) refer to this bias as the 
"finite field effect" . 

The estimation of the non-Gaussian component of C 
is non-trivial and must be derived from A^-body simu- 
lations with careful ray-tracing. Because the computa- 
tion is expensive, this covariance is commonly assumed 
to be cosmology-independent in cosmic shear studies. 
However, Eifler et al. (2009) found that covariances de- 
pend significantly on the cosmology, which impacts the 
likelihood analysis in cosmological parameter estimation. 
Currently, no ray-tracing data for such a wide range of 
cosmological parameters are available. 

Therefore, a practical method to implement this 
cosmology-dependent covariance (CDC) into one's pa- 
rameter estimation is to derive ratios between the non- 
Gaussian and Gaussian contributions from a particular 
iV-body data set and to assume that at least the ratios 
hold for other cosmological parameters. Since it is rela- 
tively inexpensive to compute the cosmology-dependent 
Gaussian covariance with Eqn. 1311 this semi-cosmology 
dependent covariance (SCDC) estimation is a useful al- 
ternative. 

Semboloni et al. (2007) first provided fitting formulae 
to enable this SCDC estimation. Sato et al. (2011; here- 
after Sll) performed a similar study but with a much 
larger data set. Hilbert et al. (2011; hereafter Hll) 
suggested log-normal approximation as a solution to the 
problem and provided detailed comparisons of their re- 
sults with those from Sato et al. (2011) and Semboloni 
et al. (2007). 

For our parameter estimation, we mainly use the fit- 
ting formulae of Sato et al. (2011) to evaluate our 
SCDC. The fitting formulae of Sato et al. (2011) do 
not provide covariance estimation for ^_ because the 
shapes are inherently complicated and cannot be accu- 
rately approximated by simple fitting formulae. The au- 
thors kindly provided a table containing these covari- 
ances for a fiducial cosmology. Consequently, in our 
likelihood analysis only C++ is cosmology-dependent. 
Cosmology-dependence is much more significant in C++ 

than C /C-i In addition, the S/N of ^+ is much 

higher than that of Therefore, modeling cosmol- 
ogy independent sample variance for ^_ (and the cross- 
covariance) is a reasonable approximation in the current 
study. This argument is supported by our experiment, 
where we perform parameter constraints with only ^+ 
data. The resulting parameter contours are highly con- 
sistent with the current case, where we use both ^+ and 

Hilbert et al. (2011) claimed that the covariance of 
Sll does not meet the positive semi-definiteness crite- 
rion, and they had to perform an eigenvalue decompo- 
sition of the Sll matrix and replace the negative eigen- 
values with zeros. However, we do not encounter this 
problem in our parameter estimation. We suspect that 
if the claim is true, perhaps our adding shot noise and 
systematic noise to the diagonal terms resolves the neg- 
ative eigenvalue problem. As the DLS consists of the 
five separate 4 sq. deg fields, we should not evaluate the 
sample variance for the ~20 sq. deg area. Instead, we 
compute the sample variance for a 4 sq. deg area and 
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Figure 24. Constraints on f7jvf and erg from DLS alone. Thie con- 
tours are obtained from ~ 20, 000 sets of parameters explored by 
our MCMC realizations. The outer and inner contours show the 
2(T and Icr boundaries, respectively. We marginalize over the un- 
certainties in the Hubble constant, photometric redshift estimation 
error, and shear calibration. 

divide the result by five. We show the matrix elements 
of the Sato et al. (2011) covariance in Figure [23l for 4 sq. 
deg area measured from the A^-body data simulated for 
the best-fit parameters of the WMAP7 study. 

4.3. Cosmological Parameter Estimation 
4.3.1. Constraints on and erg 

For the joint constraints on and erg, we only con- 
sider a flat geometry with wq ~ —1. We let ^Im and ag 
vary within the ranges [0.05,1.0] and [0.2,1.5], respec- 
tively. As mentioned above, the Hubble constant h is 
marginalized over the 0.64 < h < 0.80 range (flat prior). 
In addition, we also marginalize over the uncertainties 
both in the photometric redshift estimation ( §2.6|) and 
in the shear calibration ( §2.5.2p . 

We show the joint DLS constraints on f2jv/ and erg 
in Figure 1241 The typical "banana-shape" contours are 
seen; this is the result of the well-known degeneracy be- 
tween the two parameters in their contribution to the 
amplitude of the lensing convergence power spectrum. 
The ridge of the contours can be roughly described by 
fjg oc f^^/'^- general, contours for the Qm — fg 
constraint are highly elongated and the parameter lo- 
cation where the mimimum occurs, has little mean- 
ing. However, the relatively deep and wide DLS sur- 
vey yields a high S/N over a wide range of angle. This 
DLS cosmic shear analysis shortens the "banana" to 
an unprecedented level and enables us to constrain the 
two parameters simultaneously without relying on cosmic 
microwave background data, giving a joint constraint: 
Qm = 0.262±0.051 and fjg = 0.868±0.071. To more fully 
understand the reasons for this partial degeneracy break- 
ing, we undertake simulations using mock shear data. 

4.3.2. Impacts of Survey Depth on Cosmological Parameter 

Constraints 
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Figure 25. Survey depth vs. constraining power. We generate 
two sets of a"mock" data vector and run MCMC simula- 

tions. For a DLS-like survey, we assume Usource = 17 galaxies 
per sq. arcmin and [z) = 1.09. For the shallower survey, we use 
the parameters in the initial CFHT data studied by Hoekstra et 
al. (2006), who report Usourcfi ^ 12 galaxies per sq. arcmin and 
(z) = 0.8. The increase in the survey depth significantly improves 
the constraining power in three ways. First, the higher mean red- 
shift boosts the amplitude of the cosmic shear signal. Second, the 
higher source density reduces the shot noise. Third, the redshift- 
dependence helps to mitigate the f2jvf-o"8 degeneracy. The outer 
and inner contours show the 2cr and Icr boundaries, respectively. 

We notice that the small size of the contours for our 
joint constraints on £7 a/ and cts in Figure [24] is unprece- 
dented for cosmic shear studies. Thus, it is important 
to verify the result with simulations using mock data. 
Of course, the most robust test is to perform end-to- 
end simulations, where we create simulated DLS images 
from an input cosmology, then carry out cosmic shear 
analysis, and finally compare the results. However, be- 
cause here we are mainly interested in verifying the size 
of the contours and examining the impact of the DLS 
depth on constraining power, this full-scale simulation is 
beyond the current scope of our paper. Instead, our sim- 
ulation begins with the generation of predicted ^+ and 
^_ cosmic shear signals at the best-fit DLS parameters 
{VtM,c%)= (0.262,0.868). The main goal of the current 
simulation is to examine the impact of the DLS survey 
depth on the parameter constraining power. 

These ^+ and ^_ data points are error-free and need 
to be perturbed by the sample variance and shape noise. 
We construct a full covariance matrix from the Sll sam- 
ple variance covariance, the DLS mean galaxy number 
density, and the mean shear dispersion. We draw ran- 
dom numbers from a multivariate normal distribution 
with our full model covariance matrix. These correlated 
random numbers are added to the "exact" cosmic shear 
data vector to create the final "mock" cosmic shear data; 
the Gaussian assumption for the errors is only approxi- 
mate, and this could introduce subtle errors in the com- 
parison of the simulated data and the real data. We run 
MCMC with this mock data vector in the same way that 
we process our real DLS cosmic shear data. 

Figure [25] shows the resulting parameter constraints 



with these simulated data. The solid line is obtained 
when we assume the current depth of DLS whereas the 
dashed line results when we assume a shallower survey, 
for which we use a mean number density Usource = 12 
per sq. arcmin and the (z) = 0.8 redshift distribution re- 
ported in the first CFHT cosmic shear study (Hoekstra 
et al. 2006). This simulation verifies that the increase 
in the survey depth significantly improves the constrain- 
ing power. The reason for the S/N increase is threefold. 
First, the higher mean redshift boosts the amplitude of 
the cosmic shear signal. Second, the higher source den- 
sity reduces the shot noise. Third, the wider redshift 
baseline helps to break the JIm^cts degeneracy although 
we present here a non-tomographic study. 

A close comparison of the solid contours in Figure [25] 
with those in Figure [24] indicates a slight difference in 
size. From several runs of the above simulation, we find 
that the exact size and shape of the contour depends 
on the shape (i.e., noise from the sample variance and 
the Poissonian fluctuation) of the cosmic shear signal. It 
appears that the excess of our DLS cosmic shear signal 
at 3' < 6* < 20' mentioned in §3.2.21 may contribute to 
the overall S/N, making the contours slightly tighter. 

4.3.3. Contribution from Different Angular Scales 

Examination of Figures [T9l[2T] shows that the statistical 
errors in our cosmic shear data at > 20' are very small. 
Because our systematic errors are under good control, the 
signal on this large scale provides significant constraint 
limited only by the sample variance. To illustrate the 
point, we show two additional contours in Figure [^ 
When we use the data vectors only at 6' < 20', we cannot 
rule out the combination of an extremely high erg value 
and an extremely low ^Im value successfully. Such a com- 
bination predicts a relatively large cosmic shear signal at 
9' > 20' (this is in disagreement with our data). How- 
ever, the current shape of the DLS cosmic shear signal at 
6 < 20' favors such high as value and an extremely low 
Qm value combinations. On the other hand, if we only 
use the data at 20 < 6* < 100', it is difficult to exclude 
the low (Tg and high combination, which predicts 
a higher signal than (thus inconsistent with) the DLS 
measurements in the 9 < 20' range. Had there remained 
substantial systematics in our cosmic shear analysis, the 
contribution from the signal in the 20 < < 100' range 
would have been insignificant. 

4.3.4. Impacts of Cosmology-dependent Covariance on 
Parameter Constraints 

In the regime where cosmic shear data become sample- 
variance-limited, the role of the cosmology-dependent co- 
variance is critical. To facilitate this discussion, we show 
parameter constraints with different covariance matrices 
in Figure [27] We consider four cases here: 1) a fixed 
covariance matrix obtained by ray-tracing from the Sll 
numerical simulation [ (r2M,^^A,o'8)=(0.238, 0.762, 0.76)], 
2) a fixed covariance matrix obtained by ray-tracing 
from the Millennium Run (MR; Hilbert et al. 2009) 
[(riM,f^A,o'g)=(0.25,0.75,0.9)], 3) a fixed covariance 
matrix analytically derived by the Hll log- normal 
approximation [(17A/,f]A,crg)=(0. 279,0. 721,0.81)], 4) a 
cosmology-dependent covariance derived using the Hll 
covariance by first measuring the ratio of non- 
Gaussianity to Gaussianity and then applying it to the 
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Figure 26. Contribution of cosmic shear signals on different scales. When we use the data vectors only at 6 < 20', we cannot rule out 
the combination of an extremely high erg value and an extremely low value successfully. Such a combination predicts a relatively large 
cosmic shear signal at 9 > 20' (this is in disagreement with our data), however, without significantly violating the data points at 6 < 20'. 
On the other hand, if we only use the data at 20' < 6 < 100', it is difficult to exclude the low erg and high Qm combination, which predicts 
a higher signal than the current measurements in the 8 < 20' range. Because the data points on different scales ar c co rrelated (off-diagonal 
elements in Cs), one should not combine the results in both panels "by eye" to obtain the result shown in Figure [24l 



cosmology-dependent Gaussian covariance. For the first 
three cases where we fix the covariance matrix for a par- 
ticular cosmology, the result always favors a high-trg /low- 
J7m combination. By contrast, the SCDC implementa- 
tion with the Hll covariance in the last case gives the 
parameter constraints closely resembling the SCDC re- 
sult with the Sll covariance. As mentioned above, the 
low amplitude of the signal at 6* > 20' plays a critical role 
in disallowing low /high as combinations. However, 
Figure [27] shows that the parameter constraint can only 
benefit from this large scale cosmic shear signal when 
the covariance is allowed to vary during the MCMC run. 
One may argue that the impact of the SCDC on the 
^Ai — CTg constraint shown in Figure [27] is somewhat 
counter-intuitive because a cosmology with a large as 
might reduce and thus the SCDC might favor a large 
as cosmology. However, we find that a large as cosmol- 
ogy also substantially alters the off-diagonal elements of 
the covariance and the net effect is in fact toward in- 
creasing the value for that cosmology. 

Note that the two SCDC implementations with the 
Hll and Sll give a slight difference in contour sizes 
(lower right panel in Figure [27]) . This is because the 
covariance of Hll is higher than Sll by ~50% in the 
Gaussian regime {9 > 10') (when scaled to the same cos- 
mology) . The difference in the Gaussian regime is mainly 
due to the difference in the finite field effect correction 
between the two methods (Sll vs. Hll). Currently, no 
consensus has been reached as to which method more 
accurately corrects for the effect. 

One should not be misled into thinking that in the ab- 
sence of the cosmology-dependent covariance the result 
always favors a high-ds/low-riM combination. From our 
"mock" data simulations, we have seen opposite cases 
where the result is skewed toward a low crg/high-ilAf 
combination. For these cases, we also observe that the 
use of our SCDC tends to shift the contours in the cor- 



rect direction. However, the number of our mock data 
simulation runs is not sufficiently large to generalize this 
observation. 

4.3.5. Comparison with Other Cosmic Shear Studies 

Comparison of our DLS cosmic shear results with those 
in previous studies requires some caution. The following 
issues make it difficult to track down origins of discrepan- 
cies. First, previous cosmic shear studies obtained rela- 
tively long "banana" shapes in their VIm — o's constraints. 
Therefore, comparisons should be made at a fixed value 
of either as or fijv/ ■ Because the slope of the "banana" 
can be also different, the choice of the reference will affect 
the level of agreement. Second, theoretical improvements 
have been made in the prediction of cosmic shear signal 
and the sample variance since the first cosmic shear re- 
sults were published (Bacon et al. 2000; Wittman et al. 
2000; Van Waerbeke et al. 2000). Different methods 
of estimating the nonlinear power spectrum affect the 
normalization of the cosmic shear signal non-negligibly 
(~6% in our case). In addition, early cosmic shear stud- 
ies use Gaussian covariance matrices, which may lead to 
substantial underestimation of parameter uncertainties. 
Perhaps, counteracting this underestimation is the finite 
field effect discussed in HA.'2[ Absence of this correction 
can overestimate the parameter uncertainties. Third, dif- 
ferent authors marginalize over different ranges of nui- 
sance parameters for their target parameter estimation. 
Both the number of nuisance parameters and their ranges 
sensitively affect the parameter estimation. 

The most recent cosmic shear results are reported from 
the analysis of the SDSS Stripe 82 field by Lin et al. 
(2011) and Huff et al. (2011). The survey area is the 
largest (~275 sq. degrees) to date for cosmic shear anal- 
ysis. However, the large PSF and relatively shallow 
depth provide only moderate constraints on cosmolog- 
ical parameters. At XIm = 0.272, the result of Lin et 
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Figure 27. Impact of covariance matrix on parameter constraints. Different panels show the Qm — ag constraint results (red) for 
different covariance matrices. For easy comparison we also reproduce the result (dashed line) shown in Figure [24] which corresponds to 
the SCDC implemented using the Slf method. When we fix the covariance matrix for a particular cosmology, the result always favors a 
high-erg /low- f2ji/ combination. By contrast, if we implement the SCDC using the Hlf covariance by measuring the ratio of non-Gaussianity 
to Gaussianity and applying it to the cosmology-dependent Gaussian covariance, the resulting contours closely resemble the Sfl SCDC 
case (lower right panel). 



al. (2011) implies as = 0.63^^^. This is in excel- 
lent agreement with the resuh of Huff et ah (2011), 

who quote trg — Q-Q'^-o ib from their independent study. 
However, both results are at ~2cr tension with our value 
(erg = 0.833 ±0.034). 

Benjamin et al. (2007) performed a cosmic shear anal- 
ysis after combining the shears from the CFHTLS, RCS, 
and VIRMOS-DESCART, and GaBoDS surveys. At 
Qm = 0.24, they quote as = 0.84 ± 0.05, consistent 
with our estimate erg = 0.885 ± 0.035 for the same mat- 
ter density. From the initial 22 sq. deg CFHTLS cos- 
mic shear analysis, Hoekstra et ah (2006) estimated 
as = 0.85 ± 0.06 for flM = 0.3. Fu et al. (2007) up- 
dated the CFHTLS result using the data from a fac- 
tor of two larger area and reported CT8(f^A//0.25)*''®^ = 
0.785 ± 0.043 from their aperture-mass statistic, which 



gives as = 0.699 ± 0.038 at VLm = 0.3. Our DLS esti- 
mate as = 0.804 ± 0.021 at the same matter density is 
well within the la range of Hoekstra et al. (2006), but at 
~2a tension with the Fu et al. (2007) result. Schrabback 
et al. (2010) measured as = 0.75±0.08 at ^Im = 0.3 with 
the COSMOS tomographic analysis. This is consistent 
with our result. 

4.3.6. Joint Constraints with Cosmic Microwave 
Background Data 

Finally, we examine the cosmological parameter con- 
straints by combining the WMAP7 resultsEB ( Larson 
et al. 2011; Komatsu et al. 2011 ) and our DLS cos- 
mic shear study in a joint analysis. The contours 



1^ The MCMC data 
|http://lambda.gsfc. nasa.gov . 
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Figure 28. Joint constraints on and erg with WMAP7. The 
prior settings are identical to the ones in our DLS-only cosmological 
parameter estimation. We find as = 0.815 it 0.020 and Qa/ = 
0.278 ±0.018. For Um = 0.272 [WMAP7], we obtain o-g = 0.833 ± 
0.034, where the error is estimated with a hard prior at the quoted 
value from the DLS-only contours. 

from this combined analysis is displayed in Figure 1281 
The direction of the elongation of the contours from the 
WMAP7 results are almost perpendicular to that from 
our DLS cosmic shear results. In previous studies com- 
bining cosmic shear with the CMB data, the constraint 
from the latter is in general so strong that cosmic shear 
results in fact add only a small contribution to the joint 
constraint. However, note that the width of the "ba- 
nana" in the current study is narrower than the major 
axis of the WMAP7 ellipse. Using the same settings 
for the priors as in our DLS only experiment, we obtain 
Qm = 0.278 ± 0.018 and erg 0.815 ± 0.020 from the 
joint analysis with WMAP7. 

5. SUMMARY AND CONCLUSIONS 

We have presented a cosmic shear study from the DLS, 
which is the deepest wide-field (> 10) sq. deg cosmic 
shear survey to date. Because of the non-trivial aberra- 
tion pattern of the telescope optics, the PSF modeling 
for the stacked image was challenging. We overcome the 
difficulty by modeling the PSF exposure to exposure and 
CCD to CCD with the PCA method, and stack the PSFs 
by closely mimicking the image stacking procedure. We 
found that our initial PSF model had to be tweaked in 
two steps to remove the PSF-induced anisotropy com- 
pletely from galaxy images. The first tweak was required 
to improve the agreement between the PSF and data 
PSFs, and the second tweak removed the residual cen- 
troid bias. 

Using the star-galaxy ellipticity correlation and the B- 
mode studies, we demonstrate that the systematics in 
our cosmic shear signal are under control, which is an 
important verification before we utilize the cosmic shear 
measurement on large scales where the statistical errors 
are very small. Other systematic tests of this shape cat- 
alog using galaxy-galaxy lensing are described in Choi et 



al (2012). 

Our photometric redshifts are calibrated using > 
10,000 spectroscopic redshifts in our DLS fields. We 
employ the photometric redshift probability distribution 
p{z) instead of single-point best estimates. We show that 
the use of p(z) mitigates the systematics caused by source 
galaxies with multi-peaked, broad, or skewed photomet- 
ric redshift probability distributions. 

The shear multiplicative factor was obtained based on 
realistic image simulations. We experimented with both 
real and artificial galaxy images to derive calibration pa- 
rameters and considered the sample variance in galaxy 
population. We found that the dilution of the signal in- 
creases with magnitude, and this determines the number 
of available galaxies for lensing analysis. 

For cosmological parameter estimation, we marginal- 
ize over shear calibration error, photometric redshift un- 
certainty, and the Hubble constant. In addition, we use 
semi-cosmology-dependent covariances to avoid possible 
bias in cosmological parameter estimation when constant 
covariances at a fiducial cosmology are assumed. We ob- 
tained an unprecedentedly tight joint constraint on D,m 
and as- Without relying on other cosmological measure- 
ments, our DLS analysis alone yields Hm — 0.265 ±0.041 
and (Ts — 0.865 ± 0.077, which are consistent with the 
WMAP7 results. However, we find that when we use 
a cosmology-independent covariance matrix instead, the 
resulting parameter contours shift substantially from the 
results obtained from the use of the cosmology-dependent 
covariance. This confirms our belief that for a sample- 
variance-limited cosmic shear survey, it is critical to im- 
plement cosmology-dependence in the covariance estima- 
tion. 

Combining the current results with the WMAP7 like- 
lihood data in a joint analysis, we obtain VIm — 0.276 ± 
0.016 and as — 0.816 ± 0.019. Since our current cos- 
mic shear study is non-tomographic, we expect that our 
future tomographic study will improve this constraint. 
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